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FOREWORD 

The  work  reported  herein  describes  the  developmental  stages 
of  an  experimental  facility  to  be  used  for  investigating  the  effects 
of  injecting  a  high  current  density  into  the  dense  detonation  wave 
of  a  high  explosive.  The  experiments  represent  a  proof-of-principal 
effort  where  the  electrical  conductivity  features  of  a  detonating 
explosive  are  exploited  for  the  purpose  of  augmenting  a  number  of 
basic  explosive  properties  such  as  detonation  velocity,  detonation 
pressure  and  the  rate  of  energy  release. 

The  work  is  being  carried  out  in  the  Detonation  Physics  Branch 
of  the  Energetic  Materials  Division  of  the  Research  and  Technology 
Department  as  part  of  the  Independent  Research  Program;  Task  No. 

211,  and  NAVSEASYSCOM;  Energetic  Materials  (explosives)  Research, 
Task  No.  SRO  24-02-003/22128. 
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INTRODUCTION:  PRELIMINARY  EXPERIMENTAL  CONSIDERATIONS 

The  General  Background.  An  extensive  effort  is  underway 
to  transfer  the  stored  electrical  energy  from  a  108  kilojoule, 

(5k  V)  capacitor  bank  into  the  parti  ally- conductive  zone  behind  the 
dense  (~1C)22  particles/cm^)  detonation  front  of  an  explosive  reaction. 
The  capacitor  assembly  provides  a  versatile  and  reusable  device 
for  these  initial  feasibility  studies. 

The  principal  objective  is  to  couple  a  significant  amount 
of  this  electrical  energy  into  the  explosive  reaction  in  order  to 
effect  a  10-20%  increase  in  the  detonation  velocity.  The  subsequent 
increase  of  the  Chapin  an- Jouquet  pressure  is  expected  to  be  on  the 
order  of  20-40%. 

Figure  1  is  a  schematic  diagram  portraying  the  features 
of  the  experimental  arrangement.  Briefly,  the  108  kilojoules  of 
energy  stored  in  the  capacitor  bank  (C)  is  released  on  the  closure 
of  the  explosively-driven  closing  switch  (C.S.)  and  transferred  to 
a  storage  inductor  Ls  after  which  the  bank  is  disengaged  from 
the  circuit  via  an  opening  switch  (f).  The  use  of  the  intermediary 
storage  inductor  is  warranted  by  the  LCR  timing  constraints  associ¬ 
ated  with  coupling  into  the  high-resistive  load  typically  character¬ 
istic  of  a  detonating  explosive's  reaction  zone. 

On  referring  to  Figure  1  (a)  we  can  see  that,  during  the 
inductor-energizing  phase,  the  (unreacted)  explosive  (between  elec¬ 
trode,  E',  E"  in  the  Figure)  functions  in  the  capacity  of  a  normally- 
opened  dielectric  switch;  its  primary  purpose  is  to  hold  back  current 
from  the  discharge  circuit  (b)  until  the  inductor  has  attained  a 
maximum  current  value  of  -0.17  MA.  About  5  ysec  prior  to  the 
opening  switch  disruption  the  detonation  event  is  initiated  thus 
effectively  establishing  the  closed  LR-discharge  configuration  in 
Figure  1  (b) .  Upon  initiation  of  the  explosive,  energy  stored  in 
Ls  is  rapidly  transferred  through  the  electrodes  (~10  ys  duration) 
and  into  the  parti  ally- conductive  zone  of  the  detonation  as  it 
propagates  through  the  forward  (unreacted)  portion  of  the  solid 
explosive.  Since  energy  transmission  from  the  capacitor  bank  is 
completely  severed  from  the  system  prior  to  the  transfer  of  power 
into  the  resistive  detonation,  all  transferred  energy  is  captured 
within  the  inductive-resistive  portion  of  the  network,  and  all 
energy  is  ultimately  dissipated  within  the  resistive  load.  The 
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usual  limitations  of  impedance  matching  are  thereby  circumvented 
since  a  highly  efficient  transfer  of  energy  from  the  capacitor  bank 
to  the  storage  coil  can  be  achieved. 

The  secondary  explosive  PBX  9404  (93%  HMX)  has  initially 
been  selected  to  provide  the  source  of  the  detonation  on  which  the 
electrical  enhancement  of  detonation  velocity  (D  =  8.8  km/s)  1 
will  be  made.  These  effects  are  principally  of  an  ohmic  heating 
nature,  and  are  characterized  by  the  complex  impedance  Z  (R,wLs)  , 
the  enhancement  effect  being  tantamount  to  increasing  the  reaction 
energy  of  the  detonation  process. 

Primary  experimental  considerations  are  given  to  assessing 
the  extent  to  which  electrical  coupling  is  effected,  and  to  examining 
those  significant  transport  properties  involved  with  influencing 
such  interactions.  The  details  of  the  experimental  arrangement 
and  the  diagnostic  techniques  currently  employed  are  discussed. 

Preliminary  Experimental  Aspects.  The  experiments  represent 
a  "proof-of-principal"  attempt  at  achieving  the  transfer  of  a  large 
current  density  (on  the  order  of  ~105  amp  cm-2)  into  the  electrically 
conductive  region  immediately  forward  of  the  Chapman-Jouquet  plane 
where  the  charge  carrier  density  ne  ~  10^°  cm-2.'2  The  width  of 
this  conductive  region  (~l-3  mm)  is  controlled  by  the  chemical 
reaction  time  and  correspondingly  the  temperature.  Figure  2  depicts 
a  qualitative  profile  of  the  electrical  conductivity  distribution 
in  a  reacting  explosive.  The  more  prominent  features  such  as  the 
relative  location  of  the  conductivity  peak  and  the  extreme  exponential 
falloff  in  the  near-to-far  product  zones  have  been  experimentally 
substantiated. 2  The  spatial  distribution  indicates  the  preferred 
current  path  and  consequently  implies  information  on  the  nature 
and  extent  of  the  coupling  mechanisms  involved.  In  particular, 
the  effectiveness  of  coupling  is  far  less  significant  if  a  substan¬ 
tial  conductivity  a(e.g.,  a<_  0.1  mho  cm“l)  exist  in  the  detonation 

products.  From  a  circuit  analysis  viewpoint  where  the  conductive 
zone  is  simply  regarded  as  a  bulk  resistive  load,  the  current  transfer 
rate  is  inversely  proportional  to  the  magnitude  of  <a(£)>,  conse¬ 
quently  dictating  the  particular  network  configurations  which  are 


Dobratz,  B.  M.,  "Properties  of  Chemical  Explosives  and  Explosive 
Stimulants,"  Lawrence  Livermore  Laboratory  Report  No.  UCRL-51319, 
Revision  1,  1974. 

2 

Ershov,  A.  P.,  "Ionization  During  Detonation  of  Solid  Explosives, 
Combustion,  Explosion,  and  Shock  Waves,  pp.  798-803,  November  - 
December  1975. 

3 

Walbrecht,  E.  E.,  "Dielectric  Properties  of  Sane  Common  High 
Explosives,"  PATM  1170,  p.  1,  April  1963. 
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FIGURE  2  THE  GENERAL  FEATURES  OF  THE  ELECTRICAL  CONDUCTIVITY  EXHIBITED  IN  A 
(ONE-DIMENSIONAL)  DENSE  DETONATION  WAVE. 
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compatible  with  the  time  constraints  imposed  by  the  detonation 
velocity  through  a  given  electrode  length.  The  discussion  in  a 
later  section  will  show  that  for  a  given  explosive’s  resistance 
these  requirements  are  met  by  employing  a  simple  LR  circuit 
with  the  inductance  Ls  as  an  experimentally  adjustable  parame¬ 
ter.  The  I(t)R  drop  across  the  conduction  zone  defines  the  voltage 
requirements  imposed  on  the  storage  inductor.  In  the  course  of 
disengaging  the  capacitor  circuit  from  the  inductor,  microsecond 
switching  across  the  inductor  provides  the  necessary  voltage  required 
to  drive  the  current  through  the  resistive  load.  In  this  instance 
the  rate  determining  component  is  the  commutation  rate  of  a  paralleled 
exploding  wire  otherwise  employed  to  relax  the  voltage  stress  across 
the  opening  switch.  The  effects  of  this  current  commutation,  gener¬ 
ating  voltages  in  excess  of  102kV,  point  up  the  danger  of  overvolta¬ 
ging  various  network  components  and  in  particular  to  the  explosive 
itself.  Conversely,  it  defines  the  possible  current  limiting  condi¬ 
tions  unless  involved  modifications  to  the  electrodes  are  made. 

The  concern  over  voltage  breakdown  whether  by  channeling  through 
an  unreacted  portion  of  the  explosive  itself,  or  by  arcing  in  the 
tenuous  product  region  because  of  the  extended  conductivity  tail, 
go  into  the  considerations  of  the  particular  explosive  being  investi¬ 
gated.  For  a  given  set  of  network  conditions  it  is  necessary  that 
the  breakdown  voltage  of  the  solid  explosive  be  greater  than  any 
induced  voltages  incurred  during  the  course  of  energy  transfer  or 
premature  breakdown  at  the  explosive  may  ensue.  In  the  same  context, 
the  switching  assembly  must  be  capable  of  recovering  long  enough  to 
relax  to  the  condition  where  restrike  back  across  the  switch  is  no 
longer  a  danger. 

Both  electrical  conductivity  and  voltage  breakdown  are 
intimately  incorporated  in  the  considerations  of  the  diagnostic 
setup.  Enhancements  resulting  in  a  10-20%  increase  of  detonation 
velocity  are  shown  to  be  theoretically  feasible;4  the  actual  en¬ 
hancement  effect  will  be  measured  by  employing  fast  photography 
techniques  to  measure  velocity  changes.  The  spatial  resolution 
afforded  under  existing  streak  camera  techniques  warrants  that 
the  electrical  energy  be  roughly  comparable  to  that  of  the  chemical 
energy  released  by  the  explosive;  in  other  words,  that  some  ratio 
of  the  se  two , 


o  =  Electrical  Energy  Input 
p  "  Chemical  Release  Energy 


Unfortunately,  in  attempting  to  optimize  this  stipulation  it  is  found 
that  the  limiting  constraints  are:  1)  breakdown  voltage  requirement 
and  2)  the  coupling  interaction  time.  In  the  former  case,  there 
must  be  a  sufficient  amount  of  explosive  to  prevent  breakdown  across 
the  electrodes  during  the  course  of  switch  closure.  The  electrode 
spacing  requirement  is  even  more  stringent  when  the  concern  over 

4 

Toton,  E.  T.,  "High  Explosive  Detonation  and  Electromagnetic 
Interaction,"  NSWC  TR  79-205. 
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arcing  in  the  reaction  products  exists.  The  second  constraint 
arises  out  of  the  realization  that  the  coupling  interaction  be 
maintained  long  enough  to  provide  for  a  measurable  detonation  track. 
For  PBX  9404  with  a  detonation  velocity  of  8.8  km.  sec.“l,  it  is  found 
that  5  cm.  is  about  the  minimum  length  needed  to  establish  an  accept¬ 
able  measurement. 

The  electrodes  currently  in  use  are  designed  to  contain  a 
113.5  g,  10  cm  long  explosive  sample,  and  are  spaced  0.7  cm  apart. 

The  chemical  energy  releases  is  ~500  kilojoule.  On  the  assumption 
that  the  optimum  amount  of  electrical  energy  is  transferred,  B-0.09 
for  the  present  experimental  arrangement. 

ELECTRICAL  CONDUCTIVITY  IN  DENSE  DETONATIONS 

The  General  Significance  to  Energy  Coupling.  The  electrical 
conductivity  distribution  a(^)  plays  a  significant  role  in  determining 
(a)  the  extent  of  direct  ohmic  heating  in  the  detonation,  (b)  the 
extent  of  inductive  heating  in  the  detonation,  and  (c)  the  LR  circuit 
time-constant  characterizing  the  rate  of  energy  transfer  to  the 
detonation  wave  itself.  These  features  are  considerably  more  evident 
when  examined  in  the  context  of  the  Fourier  analyzed  LR  current 
response  where  the  first  two  Fourier  expansion  coefficients  represent 
the  dominant  energy  coupling  mechanisms.  Both  of  these  terms  are 
ultimately  responsible  for  the  ohmic  dissipation  mechanism  in  the 
conductive  region  of  the  detonation  wave,  and  both  intimately  depend 
on  the  spatial  conductivity  distribution  o  (£)  .  Details  of  the  LR 
circuit  response  are  addressed  later.  Here  a  brief  account  is  given 
of  the  physical  nature  of  a  (£ )  as  applied  to  the  coupling  mechanism 
as  a  whole.  The  displacement  relative  to  the  shock  front  is  defined 
as  £  (see  Figure  2a) . 

For  explosive  reactions  producing  a  low  carbon  content, 
the  nominal  conductivity  values  will  generally  span  a  range  of 
0.3-5  mhos  cirri. 5  Table  1  includes  the  measured  values  for  a 


5 

Brish,  A.  A.,  Tarasov,  M.  S.,  and  Tsukerman,  V.  A.,  "Electric 
Conductivity  of  the  Explosive  Products  of  Condensed  Explosives," 
Soviet  Physics  JETP,  Vol .  37,  (10),  No.  6,  pp.  1095-1099,  June  1960. 


12 


NSWC  TR  79-143 


„ 

1 

CO 

CO 

CO 

O 

0 

tu  re 

re  re 

CO 

4J 

o 

X  4-1 

X  X 

G 

re 

x  re 

mi  re 

60 

4-1  CD 

O  Q 

o  P 

G 

G 

Ml  G  G 

Pm 

Pm 

G 

C 

4-1 

o  e  x 

x 

O 

G 

G 

CO 

Ml  0 

•  s 

•  z 

U 

U 

<y 

0  o 

|3  w 

3  CP 

G 

G 

> 

44  CO 

Z 

Z 

4-4 

44 

G 

M 

a 

M 

w 

•o 

G 

od 

G 

Ptf 

> 

W 

CO 

CO 

W 

/^N 

o 

> 

B 

XI 

H 

CD  G 

PH 

CO 

X  w 

X 

H 

o 

CD  -u  CO 

w 

P 

a  a; 

o 

o 

o 

O 

Ph 

C  C  T> 

p 

CM 

CM 

r-H 

rH 

X 

d)  QJ  O 

w 

■ 

• 

• 

* 

* 

w 

4-1  G  Ml 

H 

o 

o 

o 

o 

10  u  X 

CJ 

o 

CO 

•h  4J  a 

<3 

P 

P  (U  <L) 

w 

o 

m  r-H 

M 

w 

■ 

o 

CrO 

o 

> 

{^'w'  /~N 

O  CD  CO 

■ 

+1 

o 

o 

o 

o 

pH 

■6  M  P 

o 

o 

o 

•? 

o 

^  Cd  H 

o 

r* 

o 

Cd  4-1  O 

CM 

o 

o 

CO 

a)  <-)  > 

•* 

CM 

LA 

o 

w 

Vh  O  ' 

• 

o 

A 

LA 

M 

H 

PP  > 

PQ 

H 

rH 

- 

iSH^OH 

w 

U 

Ph 

o 

g 

CM 

CM 

o 

4J 

o 

CO 

CO 

CO 

pi 

G 

•H 

PM 

bO 

4-J  Q) 

•V 

CM 

** 

j-i  cd  g 

w 

CM 

CM 

H 

rH 

CO 

LA 

LA 

LA 

LA 

X 

4-1 

Ogl-i 

> 

«-H 

r— 1 

LA 

CO 

CO 

CO 

H 

CO 

■<] 

CO 

J-f  G 

M 

. 

o 

1 — 1 

<D 

> 

o  o 

<4-1  CO 

o 

4-1 

4-1 

4-1 

44 

44 

44 

44 

44 

44 

44 

44 

44 

0 

0 

P 

G 

G 

G 

re 

re 

re 

G 

re 

G 

G 

a> 

G 

H 

w 

M 

CM 

Pi 

pi 

ps 

Pi 

oi 

0^ 

Pi 

Pi 

pi 

pi 

pi 

O 

X 

X 

Pm 

W 

/—> v 

o 

w 

G 

Pi 

v-/ 

w 

H 

> 

rH  4-1  y->. 

<c 

T3 

re 

PH 

CD 

cd  *h  6 

s 

O 

s-x 

v— ' 

O 

to 

o  >  o 

VO 

CM 

G 

•H  *H  \ 

z 

• 

O 

rH 

CO 

X 

}H  4-1  CO 

o 

o 

• 

• 

w 

QJ 

4J  G  O 

M 

rH 

rH 

P 

> 

O  3  X 

H 

•  r. 

1 

1 

CO 

J 

< 

d)  T3  0 

<3 

o 

o 

CO 

CO 

VO 

<3 

«-H  C  w 

z 

co 

lO 

i-H 

o\ 

LA 

VD 

rH 

CM 

o 

O 

o 

> 

w  o 

o 

• 

• 

• 

• 

• 

• 

1 

•H 

• 

• 

LA 

• 

u 

H 

o 

o 

o 

o 

o 

rH 

i 

o 

i-H 

n 

w 

w 

p 

CO 

<3 

fclKEM” 

CM 

CM 

00 

LA 

CM 

00 

CM 

r>- 

r^. 

LA 

rH 

< 

r^. 

i 

i 

vO 

00 

1 

1 

lO 

LA 

Ml* 

w 

mm 

« 

■  i  ■■  - 

pq 

/*s 

« 

G 

X 

s 

-5 

CO  4-ICO 

G 

v_x 

'w' 

G 

H 

V-/ 

■1 

H 

00 

V-/ 

cd  co  g 

la 

00 

vO 

»H 

00 

00 

vO 

00 

£  g  ^ 

G  W 

rH 

r-H 

o 

rH 

E 

H 

»H 

1 

1 

d 

o 

»H 

o 

Q  w 

A 

^k  ■ 

■■fill 

✓~s 

c 

44 

p 

ihB>- 

G 

cd 

s-x 

Pi 

/~s 

a) 

abw 

X 

re 

PO 

bO 

> 

53 

53 

z 

4-J 

u 

N-X 

H 

H 

H 

re 

o 

a 

X 

CO 

w 

W 

W 

B 

CM 

re 

re 

G 

O 

■9 

PM 

Ph 

P-» 

o 

o\ 

• 

bi 

SB 

S3 

nKfian 

M 

1 

+ 

O 

■*->*«« 

a 

■EP 

C 

G 

4-1 

X 

X 

H 

H 

H 

H 

LA 

in 

LA 

•H 

PQ 

G 

Z 

Z 

53 

_ a 

H 

K* 

■a 

ms 

IKS 

z 

Ph 

- 

« 

H 

H 

H 

CO 

G 

X* 

G 

O 

O 

U 

*  H 

P4 

• 

4-1 

G 

•H 

4J 

4-> 

CO 

G 

«d 

O 

G 

u 

& 

U 

G 

g 

G 

•r4 

o 

44 

CM 

u 

44 

•r4 

B 

G 

P 

G 

rG 

•r4 

4-1 

CO 

CO 

G 

44 

B 

O 

O 

o 

u 

• 

■- 

44 

(3 

>> 

s-/ 

CM 

re 

o 

bO 

i-H 

CO 

rG 

o 

rH 

4-1 

4-1 

X 

<3 

G 

•r4 

re 

s 

G 

S3 

B 

G 

G 

T) 

G 

U 

G 

o 

G 

”G 

up 

1 

co 

G 

cd 

G 

re) 

B 

G 

CO 

C 

u 

X 

re 

o 

Pm 

G 

o 

H 

G 

z 

"O 

cd 

o\ 

H 

G 

44 

1 

T3 

u 

X 

^5 

rH 

G 

« 

o 

O 

CO 

up 

S 

/■ — . 

Oi 

G 

44 

b0 

rG 

' — ^ 

in 

(3 

o 

•iM 

X 

•  iM 

T3 

13 

• 

O 

/■—s 

O 

w 

w 

X 

G 

13 

U 

re 

4-J 

X 

CO 

X 

d 

G 

u 

•H 

V4 

& 

4-1 

o 

G 

X 

G 

• 

rH 

to 

X 

G 

re 

re 

•H 

re 

T3 

P 

a 

t3 

NH 

p 

re 

G 

X 

ri 

G 

o 

G 

rH 

X 

G 

G 

PM 

O 

> 

CM 

re 

T) 

H 

G 

G 

• 

X 

G 

U 

s 

3 

G 

G 

H 

O 

cd 

co 

o 

P 

4-1 

G 

Pm 

G 

G 

•> 

G 

se 

T) 

rH 

TJ 

G 

G 

G 

U 

G 

O 

rH 

G 

Mi 

O 

G 

'G 

G 

PQ 

4-1 

G 

O 

G 

M 

<3 

P 

M 

/H 

/H 

G 

Xi 

G 

13 


NSWC  TR  79-143 


number  of  such  explosives,  but  under  rather  low  current  density 

conditions. 6*7,8, 9  Unfortunately,  the  conductivity  profile  for  PBX9404 

has  not  been  measured.*  However,  a  review  of  a  number  of  experimental 

papers  on  similar  explosives^ / 10 , 11  strongly  infers  that  the  low  current 

peak  value  will  fall  somewhere  in  the  general  range  of  those  given  in  * 

the  table.  The  mean  conductance  value  will  be  established  in  the 

course  of  our  experiment  by  monitoring  the  voltage-current  profile 

under  the  high  current  density  conditions  which  are  more  applicable  * 

to  our  particular  experiments. 

The  General  Features  of  a  Partially  Ionized  Fluid.  Ershov2 
and  Zubkov,  Luk* yanchi kov,  and  Novoselov-1-^  have  investigated  a  number 
of  general  electrical  conductivity  features  under  low  current  densities, 
i.e.,  on  the  order  of  50  A/cm2  -  specifically  for  P.E.T.N.  and 
Hexogen  (RDX)  .  Investigations  of  electrical  properties  in  explosives 
indicate  that  the  conductive  region  in  a  reactive  explosive  medium 
is  typically  characterized  by  a  3,500°K,  partially- ioni zed,  degenerate 
fluid  containing  an  exceedingly  high  density  of  neutral  molecules 
(i.e.,  no~1022  particles/cm2).2  At  these  densities  and  temperatures 
the  electrical  transport  properties  are  characteristically  attributed 
to  a  collision-dominated  mechanism.  Owing  to  the  significant  mass 
difference  between  ions  and  electrons  (me/m-[<<l)  and  to  the  pre¬ 
dominance  of  electron- neutral  collisions  (where  ve  n«1012/sec) , 
the  electrons  are  expected  to  be  the  dominant  current  carriers. 

For  the  ~3.6kV/cm  field  appearing  in  the  electrode  gap,  it  can  be 

2 

5  cm 

inferred  that  electron  current  mobilities  m  excess  10  occur. 

V  •  sec 

£ 

Jameson,  R.  L.,  Lukasik,  S.  J.,  and  Pernick,  B.  J.,  "Electrical 
Resistivity  Measurements  in  Detonating  Composition  B  and  Pentolite," 

Journal  of  Applied  Physics,  Vol.  35,  No.  3,  Part  1,  p.  714,  March  1964. 

7 

Hayes,  B. ,  "On  Electrical  Conductivity  in  Detonation  Products," 

Fourth  International  Symposium  on  Detonation,  p.  595,  1965. 

8 

Jones,  M.  S.  and  Blackman,  V.  H. ,  "Parametric  Studies  of  Explosive 
Driven  MHD  Power  Generators,"  Symposium  on  the  Engineering  Aspects 
of  MHD,  Vol.  2,  August  1964. 

9 

Jones,  M.  S.,  McKinnon,  C.  N. ,  and  Blackman,  V.  H. ,  "Generation  of 
Short  Duration  Pulses  in  Linear  MHD  Generators,"  Proceedings  of 
Fifth  Symposium  on  the  Engineering  Aspects  of  MHD,  MIT,  April  1964. 

^Zubbov,  P.  I.,  Luk 'yanchikov,  L.  A.,  and  Novoselov,  B.  S.,  "Electrical 
Conductivity  in  the  Detonation  Zone  of  Condensed  Explosives," 

Combustion,  Explosion,  and  Shock  Waves,  pp.  253-256,  June  1971. 

■^Alli son,  F.  E.,  "Detonation  Studies  in  Electric  and  Magnetic  Fields," 
Proceedings  of  the  Conference  on  Megagauss  Magnetic  Field  Generation 
by  Explosives  and  Related  Experiments,  21-23  September  1965. 

*  i 

PBX9404  (93%  HMS)  was  initially  selected  on  the  basis  of  its  very 
high  breakdown  voltage  characteristics  and  its  well  established 
detonation  properties  (cf .  Table  1)  . 
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Where  ve  n  provides  the  dominant  momentum  transfer  mechanism,  '  ' 

the  current  mobility  is  generally  expressed  in  terms  of  this  collision 
frequency  by12 


U 


e 


(2) 


where  ^  defines  the  ion-ion  collision  frequency  for  the  momentum 
transfer 'mechanism.  For  these  situations  where  the  degree  of  ioniza¬ 
tion  is  relatively  low,  the  resistivity  n(S)  is  primarily  determined 
by  electron  collisions  with  the  neutral  molecules.*  The  degree  of 
ionization,  in  turn,  depends  on  the  electron  temperature  and  the 
collision  cross-sections  involved.13  At  these  electron  densities 
the  reaction  zone  of  an  explosive  detonation  is  "too  cold"  to 
technically  be  considered  a  plasma.  Consequently,  analysis  of 
electrical  conductivity  properties,  approached  from  the  usual 
"collective-field  effect"  typically  associated  with  the  tenuous 
plasmas,  should  be  viewed  with  some  trepidation.** 

Ershov's  analysis  is  primarily  predicated  on  the  tacit 
conditions  of  low  electric  fields  and  wi dely- spaced  collisions 
(i.e.,  a  rarefield  gas)  in  the  one  limit,  and  of  low  electric  fields 
and  an  extremely  degenerate  electron  gas  (i.e.,  kT  <<3  ev)  in  the 
other. 2,15  For  dense,  weakly  ionized  fields  subjected  to  high  elec¬ 
tric  current  densities,  a  conductivity  analysis  must  generally  include 
the  "persistance-of-veloci  ty"  effects,  where  the  charge  particles 
retain  some  memory  of  past  collisions. 12  The  difficulty  ultimately 
resides  in  attempting  to  calculate  the  electron  drift  velocity 
distribution  <We>  needed  to  establish  the  electric  conductivity. 


12 

Reif,  F.,  Statistical  and  Thermal  Physics,  McGraw-Hill  Book  Company, 
New  York,  pp.  531-534,  pp.  580-582,  1965. 

13 

Dreicer,  H.,  "Electron  Collision  Phenomena  in  Partially  Ionized 
Gases,"  Phys.  Rev.  117,  p.  329,  1960. 

14 

Krall,  N.  A.,  Tnvelpiece,  A.  W. ,  Principals  of  Plasma  Physics, 
McGraw-Hill  Book  Co.,  Inc.,  New  York,  p.  57,  1973. 

15 

Zel'dovich,  Ya.  B.  and  Raizer,  Yu.  P.,  Physics  of  Shock  Waves  and 
High-Temperature  Hydrodynamic  Phenomena,  Vol.  1,  Academic  Press,  Inc., 
pp.  218-219,  pp.  229-232,  1966. 

ic 

It  can  be  shown  that  if  "self-similar"  collisions  were  frequent 
enough,  i.e.,  j  (or  ve  e)  «  ve  n,  they  would  not  effect  the 
electrical  conductivity  since  sucfi  collisions  merely  exchange  roles 

in  the  mobility  terms;  the  mobility  contributions  (where  u  e — — ) 

1  2  *aa  m  v 

a  aa 

have  been  shown  to  be  negligible  in  comparison  to  Eq.  2  above. 12,13 

*  * 

In  general,  the  plasma  "collective  field"  conditions  stipulate  that 

_3  1/2  3/2 

plasma  parameter  g<<l  where  g  =  8  x  10  (ne  '  /T '  °K) .  For  our 

detonation  process,  this  condition  is  not  satisfied  at  3500°K. 

(Even  for  ne  ~  1017/cm,  g  ~  12.). 14 
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In  the  rest  frame  of  the  detonation  the  current  mobility.  Equation 

(2),  is  related  to  the  applied  electric  field  in  the  electrode 
gap  by 

<We>=peE.  where  E  is  the  field  intensity  in  (3) 

the  gap(|E|  -  od~  s'pa5 ing  )>  where  d<*>  18  the  Potential  differ¬ 

ence.  The  corresponding  current  density  flowing  across  the  explosive's 
conductive  channel  is 


1  =  -  ne(S)e  <We>  (4) 

The  actual  mathematical  analysis  is  a  rather  involved  topic  and 
will  not  be  taken  up  here.  Dreicer's  paper1^  perhaps  best  sets  the 
stage  for  pursuing  the  quantitative  aspects  in  more  detail. 


Direct  Ohmic  Heating  Effects.  Since  it  will  be  primarily 
the  resistivity  distribution  n(£)  which  limits  the  current  flow 
through  the  detonation,  in  the  detonation  frame  of  reference,  upon 
inserting  Equation  (3)  into  Ohm's  law 

1(0  =  ne(5)  |e|yeE  ,  (5) 

and,  by  Equation  (2) 

1(C)  =  {  ne  (C )  e2/(meve,n)}E  •  (6) 

Here  ne  (£)  is  the  electron  number  density,  spatially  distributed 
throughout  the  conductive  zone,  and  defined  relative  to  the  position 
of  detonation  frpnt  £0.  Generally,  for  any  conductive  fluid  moving 
with  a  velocity  D  the  convective  current  density  is  described  by 

the  expressionJ-6 , 17 , 18 

1C  =  1  +  ne1  *  (7> 

+  It  is  readily  shown  below  that  only  the  conduction  current 

j  effects  joule  heating  term  this  increases  the  temperature 


i  a 

Jones,  M.  S.,  Bangerter,  C.  D.,  Peterson,  A.  H.,  McKinnon,  C.  D., 
"Explosive  Magneto hydrodynamics, "  AFAPL-TR  67-64,  pp.  67-72,  August 
1967. 

17 

Jackson,  J.  D.,  Classical  Electrodynamics,  2nd,  Ed.,  John  Wiley 
and  Sons,  Inc.,  pp.  471-479,  1975. 

18 

Sutton,  G.  W.  and  Sherman,  A.,  Engineering  Magnetohydrodynam ics , 
McGraw-Hill  Book  Co.,  New  York,  pp.  125-209,  and  pp.  300-304. 

1965. 
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of  the  conductive  zone  via  direct  ohmic  dissipation.  By  Equation 
(2)  and  Equation  (6)  and  n(£)  defined  by  the  expression 


n  (Z) 


(8) 


the  joule  heating  term  is  obtained.  The  energy  dissipation  in  a 
given  time  6 t  is 


Referring  to  Gauss'  Law  (or  Poisson's  equation)  gives  ^  •  1c  =  ^  • 

(ni  -  ne) .  On  the  assumption  that  E  varies  unappr eci ably  over  the 
electrode  gap  W,  the  charge  density  ne  in  Equation  (7)  scales  like 


n 


e 


e0E/W 


(10) 


The  convective  (maximum)  current  density  term  n  t)  divided  by  the 
conduction  current  density  aE  gives  ~ 


neD  *  (£oE/W)D  =  fo^  (11) 

aE  aE  “a 

For  conductivities  typically  exhibited  by  explosives  (Table  1), 
a  ~  10 2  mho/m,  W  s  0.7  x  10  2m, 
eQ  =  10  ^  farad/m,  D  =  8.8  x  102  m/sec, 

eoD  -8 

and  so  Equation  (11)  gives  — ^  -  10  .  Consequently,  for  this 

situation  the  contribution  from  the  detonation  propagation  term 
neD  to  convective  current  flow  in  Equation  (7)  is  negligible  in 
comparison  to  the  conduction  current  given  by  Equation  (6).  The 
energy  dissipation  in  a  conductive  fluid  itself  when  a  given  current 

flows  in  it,  whether  expressed  by  Equation  (9)  or  for  the  moving 
detonation  wave  (relative  to  the  lab  frame) ,  cannot  depend  on  the 
motion  of  the  conductive  system. 19  The  rate  of  evolution  of  joule 
heating  per  unit  volume  is  still  given  by  the  same  expression. 


19 

Landau,  L.  D.  and  Lxfshitz,  E.  M.,  The  Classical  Theory  of  Fields, 
2nd.  Ed.,  Addison-Wesley  Publishing  Co.,  Inc.,  Reading,  Mass., 
pp.  101-108,  1962. 
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q  =  (rj  j 2 ) 6  t .  However,  as  the  detonation  front  proceeds  through 
the  electrode  pair,  i.e.,  in  the  "laboratory  frame,"  the  measured 


"effective"  electric  field  is 


20 


E'  =  E  +  (D  x  B)/c 


(12) 


It  is  this  field  that  produces  the  conduction  current  ]  (  =  a(e)"fe') 
and  consequently  the  joule  heating  rate  per  volume  is 


,2, 


(nj*)6t  =  ( j  *E  ')  fit*  =  |jj.  (E  +  DxB)  J  fit  where  (13) 

B  is  the  magnetic  field  induced  by  current  flow  in  the  electrodes. 
If  dv  is  the  volume  of  the  conduction  zone  being  heated  by  this 
inductive  effect,  the  additional  heat  evolved  due  to  the  effect 
of  the  second  term  in  bracketed  expression  is 


q'  =  fit  j~  j  • ^DxB ^  dv 


(14) 


This  is  small  for  the  electrode  configurations  presently  being 
considered**  for  the  detonation  velocities  and  electric 
conductivities  indigenous  of  explosives.  The  details  of  this 
induced  ohmic  heating  effect  are  elaborated  on  in  the  next  Section, 
and  Appendix  A.  The  first  term  in  the  bracketed  expression.  Equation 
(13) ,  can  be  integrated  over  the  reaction  zone  volume,  then  summed 
over  each  of  the  progressive  time  increments  fit  to  give  the  amount 
of  heat  injected  into  the  detonation.  At  each  fit  time  step 


q  =  J  j.EdV  ~  ~  J  (3  •  V$  )  dV 


(15) 


j,A$  and  dV  are  all  experimentally  obtainable  quantities.  The 
potential  gradient  across  the  electrode  pair  can  be  monitored  in 
time  by  employing  standard  high  voltage  probe  techniques,  and  the 
current  flow  through  the  electrode  pair  is  measured  by  use  of  a 
modified  Rogowski  coil.  Consequently,  the  amount  of  heat  in  the 
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Landau,  L.  D.  and  Lifshitz,  E.  M. ,  Electrodynamics  of  Continuous 
Media,  Addison-Wesley  Publishing  Co.,  Inc.,  Reading,  Mass.,  p.  197- 
212,  1960. 

fhe  work  done  per  unit  time  per  unit  volume  by  the  force  field 
=  ef  E'nedv  moving  current  carriers  in  the  conduction  zone 
results  in  evolution  of  heat  and  increases  the  entropy  of  the  zone 

volume  dv  by 

r 

Some  electrode  configurations  where  inductive  heating  can  be  of  some 
significance  (e.g.,  helical  configurations)  are  readily  conceivable. 
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detonation  medium  per  unit  time,  for  each  St,  over  the  volume  of 
the  conductive  zone  results  in  q  =  <J>k  SR  (t)  for  each  time  step 

St.  With  the  subsidiary  condition  of  conservation  of  charae 
(div  j  =  0)  ,  we  can  experimentally  get 


-f 


(nj  2)dv 


and  hence  the  spatially- averaged  resistivity  for  current  can  be 
determined. 

Induced  Ohmic  Heating  Effects.  On  taking  the  curl  of 
Equation  (12),  with  E*  =  nj‘  and  inserting  Faraday's  Law  of 
Induction  for  the  (VxE)  term,  we  get  the  expression^  (reference 
is  made  to  Appendix  A  for  the  algebraic  details). 


V2B  +  V  x  (3x3),  (RMKS  units). 


The  magnetic  Reynolds  number  Rm  represents  an  indication  of  how 
effective  the  interaction  due  to  diffusion  of  the  magnetic  field 
lines  are. 

For  our  dense  detonation  situation  (Appendix  A)  , 

R  =  H0  “  (1 

m  — — 


—8  3  3 

R  ~  10  for  q  ~  10  fi*m,  D  =  8.8  x  10  m/sec,  and 

conduction  zonemlength  l  =  1.5  x  10”3  m.  t„,  the  time  it  takes 

for  the  B-field  to  dissipate  into  joule  heat,  is  given  by 


tb  =  >VL 


y  S  GC  f 


where  L  is  the  electrode  length.  With  defined  as  the  detonation 
time-of-f light  (L/D) ,  we  see  that  there  is  sufficient  time  for 
annihilation  of  the  magnetic  flux  lines  to  produce  ohmic  heating, 

i .  e .  , 


<<  1. 
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However  for  the  low  inductance  electrode  configuration,  or  (more 
significantly)  for  the  low  current  densities  which  are  available 
for  this  experiment  the  energy  release  is  about  10^  smaller  than 
that  of  the  chemical  release  energy.  Consequently,  q'  is  a  negli¬ 
gible  quantity  for  our  conditions  of  D  and  ri,  and  for  the  electrode 
configuration  which  is  currently  employed  (cf.  Equations  (A. 14), 

(A. 15) . 

BREAKDOWN  VOLTAGE  CHARACTERISTICS 


In  the  Unreacted  Explosive.  Preliminary  design  features 
of  the  experiment  had  to  address  the  problem  of  voltage  stand-off 
capability  in  various  explosives.  For  PBX  9404  (94%  HMX) ,  2  mm 
thick  samples  had  been  electrically  stressed  in  excess  of  24  kV 
without  showing  indications  of  initiation  or  voltage  breakdown 
(cf . ,  Table  1)  . 


For  nearly  instantaneous  interruption  of  the  LC  circuit 
configuration,  the  voltages  generated  are,  in  principal,  determined 
by  the  magnetic  field  diffusion  rate  across  the  storage  inductor. 

In  actual  practice  however,  the  rate  determinant  process  is  governed 
by  the  opening  switch  response  time.  To  prevent  either  the  possi¬ 
bility  of  voltage  overstressing  at  the  opening  switch  or  premature 
breakdown  at  the  explosive  (also  due  to  overvoltage) ,  an  understand¬ 
ing  of  the  particular  switch  assembly  being  employed  is  necessary. 
Also,  to  avoid  prolonged  discharge  across  the  switch  during 
the  highest  voltage  level  incurred  under  the  circuit  interruption 
impulse,  an  exploding  wire  fuse  is  employed  in  parallel  with  it; 
this  has  to  be  properly  adjusted  to  relax  the  high  voltage  stress 
generated  in  the  course  of  current  switching. 

During  the  inductor  discharge  phase,  the  peak  voltages 
across  the  explosive  coincides  with  the  peak  current;  the  induced 
voltages  generated  in  the  course  of  switching  must  satisfy  the 
restriction  where 


or  premature  breakdown  at  the  explosive  may  occur.  Here 

V^  =  the  solid  explosive's  breakdown  voltage, 

L-  =  The  equivalent  network  inductance,  including 
-1  the  storage  inductance. 

Depending  on  the  exploding  wire  length  and  material  characteristics, 
the  microsecond  switching  action  generates  voltages  in  excess  of 
100  kV. 
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In  the  Detonation  Products.  The  possibility  of  current 
channeling  through  the  product  gases  primarily  depends  on  the 
magnitude  of  the  conductivity  there.  Generally,  this  is  expected 
to  be  very  low  in  the  product  gaseous  region.  Korol'kov  et .  al. 
have  measured  high  voltage  breakdown  in  the  low  carbon  level 
detonation  products  of  PETN  and  RDX  type  explosives . 21,22  The 
values  for  voltage  breakdown  in  these  reaction  products  are  found 
to  be  comparable  to  those  generally  associated  with  solid  dielectrics. 
Korol'kov's  experimentally  observed  dielectric  strength  values  for 
RDX  are  above  65  kV  for  a  0.5  mm  gap  at  about  14-34  mm  behind  the 
shock  front.  This  would  imply  that  there  exist  nearly  a  200  kV 
holdoff  voltage  in  about  the  same  range  of  the  product  zone  for  the 
1.6  cm  electrode  gap  containing  the  PBX  9404  products  gases.*  Ref¬ 
erence  is  made  to  Figure  3.21/22  By  virtue  of  these  observed  increases 
in  dielectric  strength  during  the  1. 7-4.0  microsecond  interval 
behind  the  shock  front,  it  is  felt  that  arcing  probably  will  not 
occur.  Breakdown  in  the  far  detonation  product  zone  would  rob  the 
enhancement  process  of  the  otherwise  preferential  channeling  of 
current  through  the  reaction-to-near  product  zone  where  it  is 
most  effective.  The  object  is  to  maintain  the  voltage  across  the 
product  gases  sufficiently  low  enough  as  to  benefit  from  this  low 
conductivity  feature  of  the  product  zone.  To  further  retard  the 
tendencies  toward  arcing  in  the  product  zone,  the  copper  electrodes 
have  been  put  through  a  special  tempering  treatment  process  which 
allows  them  to  flare  outward  without  early  tearing  ahead  of  the 
reaction  products  as  the  shock  front  moves  through  them,  subsequently 
increasing  the  electrode  separation  and  lessening  chances  for  high 
voltage  breakdown.  t 

A  PRELIMINARY  CIRCUIT  ANALYSIS  OF  THE  TWO  LOOP  NETWORK 

.  Figure  1  shows  an  LCR  circuit  using  a  slow  bank  (i.e., 

I(t)~10y  A/sec)  as  the  primary  energy  source.  This  is  discharged 
into  an  intermediary  storage  inductor  which  in  turn,  serves  as  the 
means  for  the  rapid  transfer  of  energy  to  the  detonation  wave. 

Each  loop  is  examined  and  characterized  independently  with  the 
intention  of  optimizing  the  magnetic  energy  transfer  to  the  resistive 
load  shown  in  figure  1(b).  A  preliminary  examination  of  LC  storage 
circuit  (figure  1(a))  is  made  by  comparing  an  actual  current 
"ring-out"  against  a  set  of  computer  generated  derating  curves  to 
assess  the  extent  of  parasitic  components.  Analysis  is  then  made 
on  the  LR  discharge  circuit  under  the  assumption  of  rapid  switch 
commutation . 


21 

Korol'kov,  V.  L.,  Mel'nikov,  M.  A.,  and  Tsyplenko,  A.  P.,  "Dielectric 
Breakdown  of  Detonation  Products,"  Sov.  Phys.  Tech.  Phys . ,  Vol.  19, 
No.  12,  p.  1569,  June  1975. 

22 

Datsenko,  V.  A.,  Korol'kov,  V.  L. ,  Krasnov,  I.  L. ,  and  Mel'nikov, 

M.  A.,  "Dielectric  Breakdown  of  Detonation  Products,"  Sov.  Phys., 
Tech.  Phys.,  23  (1)  p.  33,  Jan  1978. 

*RDX  has  relatively  the  same  carbon  product  concentration  as  PBX  9404. 
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BREAKDOWN  VOLTAGE  FOR 
EXPLOSIVE  PRODUCTS.  (REF.  17) 
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EXPERIMENTALLY  DETERMINED  DIELECTRIC  STRENGTH  IN  THE  PRODUCT  ZONE  FOR 
PETN,  RDX  AND  90%  PETN,  AS  A  FUNCTION  OF  TIME  RELATIVE  TO  THE  SHOCK  FRONT 
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Ideal  LC  Network.  On  defining  the  ratio  of  the  surge  inpe- 
dance  to  the  loss  resistance  R  by  a  dimensionale ss  parameter  X, 
where 


lim  X  =  lim  |[  ^  £  +  00  (21) 

R*0  R+0 


is  imposed  in  the  limiting  situation  of  negligible  damping,  the 
natural  current  response  is  characterized  by  a  symmetrical,  unatten¬ 
uated  sinusoid  with  its  frequency  typically  given  by 

r  1"1/2 

a)  =  2Trf  =  <Lj  +  nLs)CJ  (Hz).  (22) 

The  inductance  in  this  instance  accounts  for  stray  component  inductances 
L.  as  well  as  for  the  storage  inductor  L  .  Specifically  if  n  defines 
the  number  of  bank  capacitors,  each  of  cfpacitance  C,  with  component 
dissipation  assumed  negligible  and  generally,  with  L  >>  L . ,  Equation 
(22)  expresses  the  characteristic  natural  frequency  <§f  the-^LC  circuit. 

V  (0)  equals  the  initial  charging  voltage  across  the  capacitor  bank. 

For  5  kV  across  the  paralleled  arrangement  of  thirty  six  240  microfarad 
(15  nH)  capacitors,  the  quarter-period  pulse  should  (ideally)  energize 
a  7.4  microhenry  inductive  load  to  a  peak  energy  of  108  kilojoules 
in  397  microseconds.  Neglecting  the  various  dissipative  losses,  the 
undamped  peak  current  condition  responsible  for  the  magnetic  flux 
energy  stored  in  the  inductor  is  simply 

1/2 

‘max  =  -Vc<nC/LS>  <23> 

=  170  kiloamperes 


In  actual  "best  case"  situations,  various  dissapative  mechanisms 
generally  reduce  this  value  by  about  20%.  The  rate  of  energy  transfer 
to  the  inductor  depends  only  on  the  initial  value  of  the  current 
derivative;  the  rate  of  current  rise  in  the  inductor  in  this  instance 
is 

It  <ni>  t-0  =  V0,/Ls  <24> 

or,  ~0.7  gigamperes/second 


LCR  Network  for  the  Capaci tor- to- Inductor  Loop.  Figure  4a 
represents  an  example  of  the  diminuted  (current)  output  from  the 
bank  "rang  down"  through  the  7.4  microhenry  inductor.  This  measured 
response  will  be  compared  against  the  damping  curves  in  Figure  4b 
to  provide  a  realistic  means  to  access  the  nature  of  the  observed 
current  suppression  accrued  by  the  circuit. 23  in  this  particular 


23 

Greenwood,  A.,  Lee,  T.  H.,  "Generalized  Damping  Cruves  and  Their 
Use  in  Solvinq  Power-Switchinq  Transients,"  Trans,  IEEE,  Vol.  82, 
Part  III,  p.  527,  1963. 
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instance  pulse  derating  is  significantly  evident  in  both  the 
amplitude  suppression  and  the  premature  quarter-period  peak; 
the  quarter-wave  misphasing  compounds  the  problem  subsequently 
resulting  in  abnormal  switching  conditions  incurred  through  switch 
mis- timing.  Once  characterization  and  adjustment  for  such  losses 
is  made,  independent  consideration  can  be  given  over  to  the  capaci¬ 
tive  switching  problems,  specifically  disconnecting  the  capacitor 
bank  from  the  storage  inductor. 


A  preliminary  analysis  of  circuit  conditions  in  the  capaci¬ 
tor-  to-  inductor  loop,  including  various  circuit  switching  components, 
indicates  that  this  more  involved  network  can  ultimately  be  reduced 
to  an  elementary  LCR  network.  The  circuits  response  is  generally 
described  by  the  second  order  equation 


or 


l  iLL  + 
L  dt  + 


i 

c 


0  , 


(25) 


(26) 


The  coefficients  Ts  and  T2  designating  the  characteristic  time 
constants  are  defined  by 


Ts  «  L/R  ,  Tp  =  RC  ; 


(27) 


inverted,  their  product 


T2  =  (L/R)  (RC)  =  LC  ,  (28) 

emerges  as  the  so-called  natural  frequency  given  by  Equation  (22) 
where 


o>o  =  (LC)"1/2  =  2irf  .  (29) 

Also,  the  quotient  expression  gives 


to  again  regain  the  condition  given  by  expression  (21)  where 
R  =  the  equivalent  network  loss  resistance, 

C  =  the  total  bank  capacity  (including  stray  capacitance), 
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L  =  the  lumped  inductances  of  the  circuit  loop 

ZQ  =  the  surge  impedence  ^  L/C  (ohms) . 

The  current  response  is  subsequently  specified  by  defining  the 
capacitor  charging  voltage  Vc(0)  and  is  generally  expressed  in  terms 
of  the  (instantaneous)  potential  driving  a  load  with  impedance  Z: 

I(t)  =  Z"1  V(t)  .  (30) 

Operational  solutions  to  Equation  (26)  are  conventionally 
expressed  by  the  inverse  Laplace  transform  response  of23 


i  (s) 


vc(°) 


s2  +  -2 —  4.  — i. 
s  T  2 
S  T 


For  the  case  where  underdamped  conditions  prevail  such  that 


R  <  2  ZQ  r  \  >  0.5  , 


we  have  that 


I(t)  =  -g 


V  (0) 


2  Ts 
2,,l/2 


(4X  “D 


(2^s)Sln( 


2_. .  1/2 


(4\  ~1) 


t/2T 


As  the  various  component  resistances  diminish  this  equation  again 
reduces  to  the  idealized  limit  of  a  simple  LC  circuit  typically 
given  by 

i  =  -V0) 

undamped  — = -  sin  w  t.  (33) 

Zo  ° 

Defining  the  dimensionless  parameter  t'  by 


t’  =  t/y[LC~r 

and  using  the  undamped  LC  case,  Equation  (33),  as  the  standard 
for  comparison.  Equation  (32)  reduces  to 


1(f) 

V^(0)/Z 


o  ( 4X  -1) 


2X 

2  . .  1/2  exp 


(^ll)  sin 


(4X 2-l)  1/2  t  *  /2X  .  (35) 
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The  family  of  curves  in  Figure  4b  are  generated  from  this.  The 
utility  of  such  curves  is  seen  in  analyzing  the  situation  depicted 
by  Figure  4a.  The  moment  the  charged  bank  is  switched,  the  ensuing 
current  is  limited  by  the  surge  impedance  of  the  network  loop  com¬ 
prising  the  bank- to- inductor  circuit  and  their  interconnections. 
Component  inductances  are  not,  a  priori,  assumed  negligible,  and 
the  inherent  circuit  resistance  further  suppresses  the  inrush  current 
response.  Observable  effects  from  stray  capacitance  (such  as  that 
between  adjacent  coil  turns  of  Ls)  are  determined  to  be  relatively 
negligible  in  this  instance.  The  most  severe  condition  is  considered 
when  the  voltage  is  at  its  peak  at  the  time  of  switching.  To  assess 
the  loop  inductance,  two  conditions  on  the  first  half-cycle  current 
waveform  are  noted.  Setting  the  first  time  derivative  of  Equation 
(32)  to  zero  gives  the  following  relation 


tan 


(4a  -1) 


1/2 


2T 


(4X2-1)1/2 


max 


(36) 


Evaluated  at  the  first  extreme  t  =  tm=x,  this  expression  gives 
the  characteristic  time  constant  Ts.  On  defining 


u  =  (4A2-1)  1/2 
“l  -  2Ts 


Equation  (36)  is  simply  recast  as 


T  =  ^  tan 
s  w. 


(w,  t  ) 
'  1  max' 


(37) 


(38) 


Also,  observing  that  the  first  crossover  point  of  the  waveform 
(t  (tt )  in  Figure  4c)  occurs  at  time 


2T 


(4A2-1)  1/2 


(39) 


makes  Equation  (38)  fully  determinable.  In  the  particular  instance 
(from  Figure  4a) 


t  (tt) 


—  =  800  microseconds 
co. 


/ 


t  =350  microseconds, 
max 
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and  we  have 


or 


(0.5)  {8xl0~4 ) 

TT 


tan 


(3 . 5x  10~4  )tt 
8xl0"4 


* 


Tg  =  640  microseconds. 


The  circuit's  loop  inductance  can  then  be  obtained  directly  by 
noting  the  onset  of  the  exponential  envelope  (i.e.,  "ultimate"  peak 
current)  is  specified  by  the  relation  (Figure  4c) 


(40) 


The  observed  first  peak  current  is  given  by 


The  realistic  network  inductance  is  then  given  in  terms  of 
(current)  waveform  trace  and  known  circuit  parameters 


(41) 


the 


% 


with 

For 


and 


L  = 


T  (TT  ) 


v°>\ 

Iobs  J 


exp 


t  (m) 

=  T  (tt)/tT  » 

vc(°) 

=  2.8  kilovolts. 

Iobs 

=  53.8  kilo amperes 

fcmax 

=  350  microseconds 

Ts 

=  640  microseconds 

(42) 


The  circuit  loop  inductance  is  17.4  microhenries;  subse¬ 
quently  stray  and  component  inductances  are  particularly  significant 
in  this  case.  Equation  (23)  for  peak  current  with  no  damping  gives 


2,800  V-<,36H,240) 


=  62,392  amperes  peak 
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To  effect  the  peak  current  portrayed  by  the  current  waveform  (in 
Figure  4a),  a  per  unit  current  reduction  of  53.8/62.4  =  .86  must 
occur.  Referring  to  Figure  4b  indicates  that  X  =  7  satisfies  our 
conditions.  With  the  surge  inpedance  Z  =  45  milliohms,  the  loop 
loss  resistance  as  given  by 

R  =  Z  /X 

(43) 

-  6.  4  m Q 

This  is  the  resistance  limiting  the  current  to  53,800  amperes  in 
the  capaci tor-to-inductor  loop  otherwise  approaching  55,748  amperes. 

The  instantaneous  power  transferred  is  resolved  into  active  (i.e., 
dissapative  or  energy  absorbing)  and  reactive  (i.e.,  energy-storing) 
components;  the  respective  phase  factors  are  the  power  factor  and 
the  reactive  factor. 


Addressing  for  the  moment  only  the  dissapative  component,  the  time 
average  power  dissipated  to  the  circuit  during  the  first  quarter- 
period  is  typically 


<P> 


1 


t 


max 


0 


t 


max 


I(t)Vc(t)dt  . 


(44) 


From  this  the  "effective  power"  is  given  by 
I  I  it) 


P  = 


max 


max 


Vc  (0)  x  (p.  f.) 


(45) 


The  extent  of  power  dissipation  in  the  circuit  during  the  inductor 
energizing  phase  can  be  illustrated  by  the  power  factor  (p.f. ).  A 
low  dissipation  factor  is  a  good  criterion  for  the  circuit  perform¬ 
ance  here.  For  a  series,  inductor  dominated,  LCR  circuit  the  (lag) 
power  factor  is  given  by 


(P.f.)  = 


Tfr 


+  (to^L 


-  1  ) 


(46) 


(the  active  component  phase  factor  between 
the  instanteous  current  and  the  instantaneous 
voltage) 


where  L  =  L  +  S L  . 

s  3  3 
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As  (p.f.)  decreases,  progressively  greater  currents  are  delivered  to 
the  inductor  L  .  The  ratio  of  the  maximum  energy  stored  to  the 
energy  loss  ina  given  power  pulse,  conventionally  referred  to  as 
the  quality  factor  (Q)  for  the  circuit,  is  a  directly  measurable 
quantity  that  indicates  the  extent  of  damping  and  time  misphasing. 
For  a  series  LCR  circuit  this  is  generally  defined  as 


Q  =  0)oL/R 


(47) 


Consequently,  the  reactive  factor  Q  (or  equivalently,  X)  can  be  used 
to  explicitly  determine  the  amount  of  stray  inductance  L.,  and  hence 
accurately  determines  the  extent  of  misphasing  in  the  circuit.  With 


Figure  4d  illustrates  the  generalized  family  of  power  factor 
curves  that  one  can  typically  generate  directly  from  equation  (49) 
with  the  characteristic  parameter  being  Q.  As  a  rather  dramatic 
example  of  power  lost,  initially  where  53.8  kiloamperes  is  trans¬ 
ferred  under  "lossy"  conditions:  ,X  =  7;  p.f.  =  0.4.  From  Equation 
(45)  the  "effective  power"  pulse  dissipated  is  .17  gigawatts  or  about 
40%  of  the  total  available  power  pulse. 


For  systems  with  large  Q  and  X,  u.  =  o)Q  as  Figures  4b  and  4d 
indicate  the  time  misphasing  diminishes  andef f ect i ve  power  comes 
principally  from  the  reactive  (i.e.,  the  energy  storage)  component 
of  instantaneous  power.  The  optimum  condition  for  the  LR  circuit 
is  to  have  a  high  Q  and  a  low  power  factor.  This  optimum  condition 
can  be  approached  in  the  present  facility  by  using  circuit  lines 
with  the  lowest  practical  resistance  and  length. 


LR  Discharge  Network.  Current-detonation  interaction  commences 
with  the  disengagement  of  the  capacitor  bank  and  switching  assembly. 
About  5  microseconds  prior  to  the  opening  switch  disruption  the  HE 
detonation  event  is  initiated,  thus  effectively  establishing  the 
closed  LR-discharge  configuration.  From  this  point  the  conductive 
channel  of  the  reacting  explosive  serves  to  complete  the  current 
path  across  the  explosive-electrode  pair.  By  a  proper  adjustment 
of  the  storage  inductor  to  match  the  characteristic  resistance  of 
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this  conductive  channel  the  resulting  ratio  of  (R/Le)  uniquely 
prescribes  the  decay  rate  of  the  inductively  stored1" energy  in  Lg  and 
the  subsequent  rate  of  energy  dissipation  to  the  detonation  wave. 


The  current  response  from  a  series  LR  discharge  circuit 
is  typically  described  by  the  circuit  impedance  Z  =  R  -  iwLs. 

Figure  5  illustrates  the  current  transfer  to  the  conduct ionzone 
of  the  detonation  wave  for  a  set  of  detonation  conductivities  ranging 
from  0.4  to  4.0  mho/cm.  L  is  the  experimentally  adjustable  paramete 
required  to  match  a  given  resistance.  At  this  point  of  events  the 
capacitor  bank  has  been  decoupled  from  the  network  (see  Figure  la,  b) 
The  case  illustrated  in  Figure  5  uses  a  storage  inductive  value  of 
7.4  microhenries.  Any  decrease  of  magnetic  field  energy  from  the 
storage  inductor  must  ultimately  show  up  as  dissipated  energy  in  the 
resistive  components  of  the  circuit.  Consequently  the  energy  pre¬ 
dominately  goes  into  the  conduction  zone  of  the  detonation  wave  as 
joule  heat,  i.e.,  n(C)j^  (c.f.  equation  8).  On  commencement  of  the 
detonation  (designated  as  a  time  td)  the  current  response  is 
described  by  conservation  of  energy  as 


d_ 

dt 


L  I2  (K,  t) 

— - ^ -  =  1 (Cft)R(C),  for  t  >  td. 


(50) 


or 


Ls  +  IR  =  0. 


(51) 


with  the  definition  y= 

N 

linear  homogeneous  D.E. 


R(£) 

L_ 


Equation  (51)  reduces  to  the  simple 


k  + 


KCft)  =  0, 


where  the  explosive's  resistance  is  R(£). 


(52) 


Explicitly  noting  that  the  current  transferred  through  the  conduction 
zone  is  position  dependent  reiterates  the  fact  conveyed  earlier  that 
the  conduction  zone  is  a  nonisotropic,  resistive  fluid  with  exponen¬ 
tially  increasing  resistivity  in  the  product  region.  This  position 
dependence  of  conductivity  is  not  used  in  the  general  treatment 
presented  in  this  report.  However,  to  treat  the  details  of  energy 
released  in  the  reaction  zone  and  how  this  energy  is  communicated  to 
the  detonation  front  will  require  knowing  this  spatial  relationship 
of  conduct i vi ty . 
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FIGURE  5  THE  IDEALIZED  CURRENT  TRANSFERRED  FROM  THE  STORAGE  INDUCTOR  TO  THE 
CONDUCTION  ZONE  OF  THE  DETONATION  WAVE.  FOR  VARIOUS  EXPLOSIVE 
CONDUCTIVITIES,  RANGING  FROM  0.4  TO  4.0  mho/cm.,  L$  IS  THE  EXPERIMENTALLY 
ADJUSTABLE  PARAMETER.  THE  CAPACITOR  BANK  HAS  BEEN  DECOUPLED  FROM 
THE  NETWORK. 
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Equation  (52)  can,  of  course,  be  directly  integrated,  and 
on  employing  the  initial  condition  I(xD)  =  I  we  get  the  familiar 
exponential  decay  of  I  in  the  storage  inductor  with  time: 

I(£,t)  =  I0  e  for  t>xD,  with  the  decrement  (53) 

being 

y  =  if  L  >>  Z  L.  .  (54) 

*  s  .  i 

a  n  J 


The  natural  current  response  is  specified  by  Equation  (53)  for  times 
t  >>  t  and  implies  inf  ini  tely-fast  switching.  In  actual  practice, 
however,  the  current  is  driven  by  the  (~2  ys)  voltage  spike  generated 
during  fuse  rupture.  The  forcing  function  superimposes  an  additional 
"particular"  solution.  Consequently,  in  the  time  domain,  the  switch¬ 
ing  surge  produces  a  "forced"  response  which  is  given  by 

/'-Rt(Lj+Ls)/L 


I(t)  = 


-Vo 


1 


-exp 


V 


R 


1+  tVv^vtLj+v]* 


exp  (-yt) 


(55) 


Under  the  assumption  that  the  fuse-assembly  inductance,  L.  <<  L  , 
Equation  (55)  reduces  to  the  abbreviated  form  given  by  Eqelations  (53)  . 
The  particular  scheme  used  here  thus  selects  out  the  peak  current 
I  (td)  =  I  excited  at  the  frequency  y  which  represents  the  natural 
response  If  L  is  properly  adjusted  to  the  resistance  of  a  specific 
explosive.  Tnis  establishes  the  value  for  y  which  corresponds  to  the 
transit  time  through  the  electrodes  -  the  latter  being  preselected  on 
the  basis  of  the  3-criteria  established  earlier  by  Equation  (1.1). 

The  current  response  for  a  typical  range  of  (low  current)  explosive 
resistances  in  the  time  domain  Equation  (53)  is  shown  in  Figure  5. 


Dissipative  Power  Available  to  the  Explosive.  Under  actual 
power  transfer  conditions,  the  power  absorbed  by  the  modular  opening 
switch  unit  itself  must  be  accounted  for.  Some  additional  loss  to 
vaporization  of  the  exploding  fuze  (see  circuit  in  Figure  5)  also 
depletes  a  fraction  of  the  available  commutated  energy  at  a  rate 


WQ  =  d_  (LI2)  .  (56) 

dt  2 

The  remaining  power  available  to  the  explosive  detonation  wave,  on 
being  driven  by  the  circuit  breaking  unit's  voltage  is 

W  =  IV  - 

-  10  gigawatts. 
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Wo  is  the  power  loss  to  the  circuit  breaking  unit. 

.  If  stray  inductance  effects  and  switch  losses  are 
assumed  relatively  negligible  compared  to  L  and  the  detonation  wave 
resistance,  the  voltage  for  the  LR-discharge  circuit  being  driven 
by  the  (parallel)  switching  unit  is  inferred  (under  quasi-static 

conditions^)  as 


V  =  IZ'  (59) 

where  the  real  part  of  Z  is  Z1  and  the  imaginary  part  is  Z'',  i.e.. 


Z  =  Z'  +  iZ"  (60) 

The  expression  for  Z  can  be  regarded  as  representing  the  first  terms 
in  an  expansion  of  the  functions  Z(y)  in  powers  of  the  frequency. 

The  first  term  Z*  represents  a  pure  resistance,  hence  direct  ohmic 
heating.  The  effect  due  to  the  second  term  Z'*,  arises  because 
the  magnetic  field  induced  by  a  variable  current  in  the  electrodes 
generates  eddy  currents  and  therefore  an  additional  dissipation 
of  energy  in  the  detonation. 20 , 24  This  "induced"  ohmic  effect  is 
discussed  on  page  17  from  a  different  viewpoint  and  detailed  in 
Appendix  A. 

The  joule  heating  (i.e.,  the  power  dissipated)  in  the  detonation  wave 
conduction  zone  is  purely  D.C.  ohmic  in  nature  and  is  thus  simply 
given  by 

I2Z*  =  I2R(U 

where  R(£)  is  the  conduction  zone  resistance.  I  is  the  current 
transferred  through  that  resistance  when  driven  by  the  non-linear 
EMF  generated  by  the  exploding  wire  fuse.  Within  the  particular 
idealized  experimental  limitations  considered  here  (eg.  no  switch 
losses,  etc.)  the  maximum  rate  of  energy  transfer  from  the  7.4  yH 
storage  inductor  is  V/L  =  13.5  GW  of  power  driven  into  the  explosive 
load.  s 


For  a  L  /R  of  0.4ysec  (See  Figure  5)  and  a  postulated  trans¬ 
fer  of  170  kA,  the  maximum  energy  transfer  to  the  explosive's  detonation 
wave  is  found  by  integrating  the  instantaneous  power  expression  over 
time.  This  results  in 


[  exp  (-2  2-  TD) 

=  107  kilojoules 


(61) 


24 

Chu,  B.  T. ,  "Thermodynamics  of  Electrically  Conducting  Fluids," 
Physics  of  Fluids,  Vol.  2,  No.  5,  p.  473,  Oct  1959. 
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For  a  10  gm  explosive  sample  with  approximately  5  KJ/g  detonation 
energy  density,  there  could  be  as  much  as  a  factor  of  2.2  increase 
in  observed  enhancement  of  energy  for  the  present  experimental 
arrangement. 


As  calculated  in  Appendix  A,  for  the  particular 
coaxial  accelerator  geometry  and  explosive  resistivities  considered 
here,  the  101  J  incurred  through  magnetic  diffusion  is  insignificant 
compared  to  the  direct-contact  ohmic  heating  effect  of  Equation  61. 
Under  properly  parameterized  conditions  the  contribution  to  deposited 
energy  by  the  magnetic  diffusion  process  can  be  large.  This  subject 
will  be  addressed  in  future  work. 

A  DESCRIPTION  OF  THE  EXPERIMENTAL  ARRANGEMENT  USING  A  FIXED 
CAPACITIVE  ENERGY 


The  actual  experimental  facility  is  shown  in  Figure  6.  It 
was  schematically  depicted  in  Figure  1  as  a  two  loop  network  com¬ 
prised  of  the  bank  array,  the  storage  inductor,  the  resistive  load 
(i.e.,  the  reacting  explosive),  and  various  switching  assemblies. 

In  actual  practice  these  circuits  entail  a  number  of  surge  components 
comprising  a  variety  of  dissipative  and  transient  effects.  A  scheme 
employing  a  set  of  generalized  damping  curves  is  used  to  characterize 
the  LC  loop  (shown  in  Figure  la)  under  actual  network  loading  condi¬ 
tions  . 


The  capacitor  bank  currently  in  use  consists  of  nine 
modules,  each  containing  four,  240  microfarad  (5  kV)  capacitor  units. 
The  bank  array  is  configured  in  a  parallel  arrangement  to  deliver 
optimal  current  to  a  7.4  microhenry  inductive  load  during  the  first 
quarter  period  of  a  630  hertz,  underdamped  osci  11a  t’*  .  Ideally, 

the  waveform  is  truncated  at  the  first  current  maximum  to  achieve 
optimal  energy  density  at  the  inductor.  This  is  accomplished  by 
rapidly  interrupting  the  capaci tor- to- inductor  circuit  at  that 
instant  -  the  generated  emf  across  the  inductor  going  into  driving 
the  current  through  the  resistive  (explosive)  load  hence  transferring 
the  stored  magnetic  energy  to  the  detonation  wave. 

Direct  coupling  from  the  slow  bank  to  a  highly  resistive 
explosive  produces  an  unfavorably  long  pulse,  subsequently  imparting 
only  a  small  fraction  of  the  available  capacitive  energy  into  the 
detonation.  The  nominal  resistive  load  values  typically  character¬ 
istic  of  the  dense  detonation  processes  envisioned  here  indicate 
that  it  is  not  possible  to  increase  the  RC-time  constant  much  above 
a  few  milliseconds.  However,  satisfying  Equation  (1)  requires  that 
the  detonation  front  transit  through  the  electrode  pair  in  about 
10  microseconds  -  the  current  pulse  thus  lags  by  about  two  orders  of 
magnitude.  To  circumvent  this  problem  the  low  impedence  storage 
inductor  is  inserted  in  parallel  between  the  bank  and  the  electrode 
pair  to  serve  as  a  fast  response,  intermediary  current  source. 

The  Use  of  Inductive  Storage  Schemes  for  Rapid  Energy  Trans- 
fer.  Salqe  and  Braunsberger  have  discussed  the  use  of  various 
inductive  storage  schemes  to  pulse  shape  the  current  output  from 
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FIGURE  6  ACTUAL  TEST  FACILITY  SETUP 
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large  capacitor  hanks,  specifically  in  applications  related  to  high 
amplitude,  microsecond  current  responses. 25  Once  the  mean,  high 
current  resistance  of  the  detonation  wave  has  been  established  an 
adjustable,  bitter-type  storage  inductor  can  be  employed  to  replace 
the  existing  7.4  microhenry  helical  unit.  This  will  provide  the 
necessary  versatility  to  tune  the  network  inductance  for  optimizing 
the  capacitor  energy  transfer  to  the  resistive  detonation  wave. 

During  the  course  of  designing  this  experiment,  two  helical 
solenoid  designs  were  employed,  both  using  high  quality  copper  stock 
as  the  conductive  element.  The  conversion  of  a  high  energy  capacitor 
pulse  of  0.17  MA  generates  an  axial  magnetic  field  near  a  quarter 
of  a  megagauss.  Consequently,  mechanical  stress  loads  on  the  order 
of  246  MPa  are  realizable.  The  initial  design  given  in  Figure  7 
consist  of  a  0.5  inch  solid  stock,  10-turn  helical  coil,  internally 
constrained  by  a  4.5  inch  diameter  "Lexan"  mandrel  to  check  the 
inward  radial  over-stressing,  and  the  compressive  forces  of  the  end 
turns.  A  highly  dielectric  elastomer  separated  the  individual 
windings.  The  entire  coil  assembly  was  encapsulated  by  a  one  inch 
thick,  epoxy  -  reinforced,  fiber-glass  boot.  This  construction  is 
very  bulky  and  lacks  required  versatility  when  changing  the  inductance 
value,  i.e.,  by  installing  additional  turns,  or  coils.  The  second 
design,  the  Bitter  Disks  type  construction  is  also  a  single  layer 
solenoid.  However,  instead  of  filamentary  type  windings,  this 
design  employs  a  set  of  copper  disks;  each  turn  is  split  radially  and 
the  split  edge  slightly  displaced  out  of  the  disk  planeto  form  a 
step  which  overlaps  with  the  previous,  adjacent  disk.  "The  disks  are 
interleaved  with  each  other  and  with  similar  disks  of  insulation  to 
form  a  continuous  helix",  as  illustrated  in  Figure  8.26  Clamping 
the  desired  number  of  disk  stacks  together  provides  a  versatile  and 
relatively  compact  storage  inductor. 

Component  and  Switching  Assemblies. 

ARRANGEMENT  FOR  SEQUENCING  EVENTS.  Figure  9  is  a 
block  diagram  illustrating  the  various  events  as  they  typically  occur 
in  the  actual  experimental  setup.  Trigger  information  via  a  central 
control  console  (ITIG)  enjoins  three  main-channel  networks  to  properly 
coordinate  the  respective  timing  of  the  various  firing  line  and 
switching  assemblies.  As  depicted  in  Figure  9  (by  the  dashed  network) 
a  357  microsecond  delayed  trigger  pulse  to  the  Explosive  Initiation 
Network  (EIN)  evolves  from  this  point.  EIN  enlists  the  use  of  two 
10  kV  pulsers,  one  of  which  initiates  an  exploding  bridgewire  deton¬ 
ator  (E.  B.  W.  DETR)  which  in  turn  initiates  a  3  cm  pentolite  booster 

25 

Salge,  J. ,  and  Braunsberger ,  U.,  "Inductive  Energy  Storage  Systems 
Applied  for  the  Extension  of  Current  Pulse-Duration  of  Capacitor 
Banks,"  Symp.  Fusion  Technology,  7th,  Grenoble,  Oct  1972. 

26 

Frungel,  F.  B.  A.,  High  Speed  Pulse  Technology,  Vol.  3,  Academic 
Press,  Inc.,  New  York,  p.  353,  1976. 
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ASSEMBLY  COMPONENTS  OF  THE  1.0  AND  7.4pH  HELICAL  STORAGE  INDUCTORS. 
Clockwise;  7.4  /iH  coil  with  "Lexan"  Mandrel  inserted;  coil  system  shown  encapsulated  in  a  rigid, 
high-dielectric  elastomer;  assembly  inserted  in  a  tight  fitting  1.5  inch  thick,  high  tensile  strength, 
fiberglass-epoxy  composite  shoe.  Below  the  7.4  ^H  assembly  is  a  1.0  jttH,  multi-turn  coil 
unconstrained. 


E  F 

FIGURE  8  SCHEMATIC  DIAGRAM  DEPICTING  COMPONENTS  FOR 
A  BITTER-TYPE,  DISK  INDUCTOR 
Coil  Components: 

(A)  stamped  conductor  disk  with  insulator  disk  (B) 

(C)  milled  conductor  disk  with  insulator  disk  (D) 

(E)  assembly,  side  view 

(F)  assembly,  cross  section 
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interfacing  the  HE  sample.  The  other  unit  307  microseconds  delayed, 
supplies  a  103  7  kv  pulse  to  a  3  cm  (2  mil  tungsten)  exploding  bridge- 
wire  (EBW)  light  source.  The  100  microsecond  light  pulse  generated 
from  this  EBW  provides  an  intense  frontal  illumination  source  for 
streak  camera  photography.  The  EIN  is  initially  delayed  307  micro¬ 
seconds  to  allow  adequate  time  for  the  electrical  response  of  the 
opening  switch.  Terminating  the  EIN  is  113  grams  of  explosive. 

The  second  information  channel  initiates  the  Si-trigger 
circuit.  activates  the  switching  assembly  to  clamp  the  capacitor 

bank  on  line  once  it  is  up  to  charge.  A  7  kV  transformer  isolated 
pulser  initiates  three  (EBW)  detonators  to  explosively  drive  the 
clamp  switch  into  closure.  The  enjoining  circuit  (in  Figure  5) 
is  further  delayed  317  microseconds  to  allow  for  the  current  buildup 
in  the  storage  inductor.  Time-based  from  the  clamp  switch  closure 
the  trigger  pulse  S2  is  retarded  90  microseconds  to  allow  the 
modular  opening  switch  (MOS)  the  necessary  induction  time  required 
for  full  mechanical  response  and  a  subsequent  40  microseconds  to 
complete  voltage  recovery.  Figure  10  depicts  the  idealized  electrical 
responses  in  various  circuits  displayed  in  the  correct  timing  sequence. 

27  28 

DESCRIPTION  OF  THE  SWITCH  ASSEMBLIES.  '  The 
switch  assemblies  currently  employed  are  part  of  a  developed  tech¬ 
nology  used  in  the  Naval  Research  Laboratory's  "Linus"  fusion  reactor 
program.  Only  a  brief  description  of  these  are  presented  here.  The 
capacitor  bank  is  brought  on  line  by  rupturing  a  solid-dielectric 
insulator  such  as  that  employed  by  the  switching  arrangement  illus¬ 
trated  in  Figure  11.  The  shutter  mechanism  for  this  switching  unit 
enlists  the  use  of  a  217  mg  explosive  EBW  detonator  to  mechanically 
drive  a  small  section  of  aluminum  channel  plate  through  the  dielectric 
stand-off  (as  illustrated  in  Figure  11a) .  The  plate  is  connected  to 
the  low  side  of  the  bank  circuit.  Experiments  using  various  thick¬ 
nesses  of  polyethylene  indicate  that  the  driven  section  channel  (a 
2-3  cm  cylindrical  boss)  is  capable  of  completing  closure  with  about 
+  0.5  microsecond  jitter.  To  deliver  the  (5  kV)  170  kA  current, 
three  switch  units  are  fired  in  parallel.  On  a  trigger  command 
these  are  simultaneously  driven  through  the  polyethylene  dielectric 
and  into  a  set  of  contact  anvils  (Figure  lib).  The  advantages  of 
these  solid-dielectric  switches  over  other  types  of  conventional 
closure  devices  are  derived  from  the  inherent  simplicity  of  a  unit 
with  the  capability  of  closing  a  circuit  with  a  large,  slow  pulse 
and  yet  maintaining  minimal  shutter  jitter,  pulse  distortion  and 
negligible  power  loss.  Using  a  0.05  cm  thick  piece  of  polyethylene 
to  hold  off  5  kV  the  observed  response  time  for  three  paralleled 
switch  units  is  about  12.0  +  0.5  microseconds. 


27 

Ford,  R.  D.,  and  Young,  M.  P.,  "Chemical  Detonator  Solid-Dielectric 
Switches  for  Starting  and  Clamp  Applications,"  Proceedings  of  the 
8th  Symp.  on  Fusion  Tech.,  EUR  5182e,  Sep  1974. 

28 

Ford,  R.  D.,  and  Vithovitsky,  I.  M. ,  "Explosively  Activated  100  kA 
Opening  Switch  for  High  Voltage  Application,"  NRL  Report  No.  3561, 
Jul  1977. 
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FIGURE  10  IDEALIZED  SWITCH  RESPONSE  PROFILES:  TIMING  SEQUENCE  OF  EVENTS 
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FIGURE  11  SOLID-DIELECTRIC  CLAMP  SWITCH  USED  IN  CAPACITOR  STARTING 
APPLICATION.  (REF.  27) 
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The  bank  is  disengaged  from  the  circuit  by  mechanically 
severing  the  circuit  loop.  The  opening  switch  shown  in  Figure  12 
is  a  modular  assembly  designed  to  interrupt  the  current  in  the 
capicator-to-inductor  loop  at  peak  current.  The  particular  breaker 
switch  presently  being  used  is  an  explosively  activated,  twelve 
stage  switch  capable  of  achieving  the  necessary  electrical  inter¬ 
ruption  in  less  than  40  microseconds.  Parallelling  a  60  cm  fila¬ 
mentary  exploding  wire  in  a  low  inductance  configuration  with  the 
switch  provides  the  required  current  commutation  rate  for  necessary 
voltage  recover,  consequently  suppressing  arcs  and  switch  restrike 
to  individual  switch  sections.  Under  optimal  conditions  the  power 
consumption  required  to  vaporize  these  copper  wire  fuses  is  expected 
to  be  relatively  low  -  depending  only  on  the  mass  and  heat  of 
vaporization  of  the  copper.  The  use  of  aluminum  foil  fuses  in  a 
water  tampered  arrangement  results  in  lower  power  expenditure. 

The  water  medium  also  offers  the  advantage  of  generating  up  to 
three  times  the  driving  voltage  compared  to  that  in  air.  The 
dictating  parameters  for  power  consumption  here  are  the  rate  of 
recovery  of  the  switch  and  the  amplitude  of  the  commutated  current 
during  commutation  in  the  fuse. 

In  the  closed  circuit  posture  the  conducting  switch 
element  is  a  12  inch  long,  cylindrical  aluminum  liner.  The  liner 
is  completely  filled  with  a  2  inch  diameter  solid  paraffin  plug  - 
along  the  entire  axial  length  of  which  runs  a  150  grain  piece  of 
PETN  "pr ima-cord. "  The  cartridge  is  connected  in  series  to  the 
low  voltage  side  of  the  clamp  switch  on  the  one  end  and  the  high 
side  of  the  inductive  store  on  the  other.  On  detonation  of  the 
pr  ima-cord,  the  action  of  hot  explosive  gases  is  transferred  to 
the  paraffin.  During  the  opening  process  the  high  speed  molten 
paraffin  issues  between  a  set  of  tapered  polyethylene  insulating 
rings  to  quench  the  high  voltage  arcing.  After  opening  over  the 
bending  dies  each  stage  of  the  switch  is  separated  from  another  by 
the  polyethylene  rings;  these  partitions  serve  to  insulate  the 
respective  sections  of  the  segmented  liner.  The  qualitative  profile 
shown  in  Figure  10  portrays  a  mechanical  response  time  (i.e.,  the 
explosive  induction  time,  plus  a  time  for  elastic  yielding  of  the 
aluminum  liner  over  the  bending  dies)  of  90  microseconds.  The  arc 
dwell  .time  for  completing  current  interruption  depends  on  the  cross 
sectional  area  and  inductive  behavior  of  the  particular  fuse  con¬ 
figuration  paralleled  across  the  switch. 

Other  conventional  schemes  generally  employed  as  high  power 
circuit  breakers,  such  as  ignitron  actuated  "power -crowbar"  techniques 
or  ohmically  vaporized  shunts  generally  require  a  rather  substantial 
power  consumption.  Both  of  these  techniques  are  believed  to  be 
inadequate  for  the  current  parameters  considered  here.  Ignitron 
circuits  are  generally  rather  expensive,  electronically  complex  arrays 
which  can  induce  significant  spurious  noise  into  the  current  profile. 
Vaporizing  shunt  devices  can  expend  upwards  of  50%  of  the  system's 
initially  carry  the  large  currents. 
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ELECTRODE  ACCELERATOR  CONFIGURATIONS.  Figure  13  illustrates 
the  two  electrode  configurations  currently  considered  for  this  study. 
In  each  instance  the  explosive  sample  is  extended  2.5  cm  beyond  the 
electrode  ends  to  insure  that  a  planar  detonation  front  is  propagating 
as  nearly  normal  to  incident  to  the  electrode  entry  as  is  possible. 

The  explosive  extensions  also  serve  to  establish  baselines  for  the 
camera  measurements.  The  planar  configuration  affords  a  simpler 
analysis  because  the  electric  field  distribution  in  the  rectangular 
shaped  reaction  zone  is  uniform  and  one  dimensional.  For  this  con¬ 
figuration  the  camera  will  be  able  to  directly  track  the  detonation 
front  down  the  length  of  explosive  slab  using  standard  optical  tech¬ 
niques.  The  principal  concern  with  this  geometry  lies  with  the  risk 
of  "current  channeling"  (arcing)  along  the  boundary  walls  of  the 
accelerator.  The  electrode  configuration  to  be  employed  first  is  a 
coaxial  geometry.  Here  the  explosive  conduction  path  is  an  annular 
channel.  The  obvious  disadvantage  here  lies  in  the  complication  of 
the  electric  field  distribution  which  is  radially  divergent.  However, 
this  geometry  lends  itself  to  direct  velocity  measurements.  (These 
are  of  a  similar  nature  to  established  cylinder  tests.)  The  annular 
geometry  (Figure  14)  does  not  present  conditions  for  current  arcing  to 
occur.  To  assess  detonation  velocity  enhancement,  a  high  speed 
image-converter  camera  is  employed  to  track  the  expansion  rate  of  the 
outer  cylinder  wall.  These  results  are  compared  against  the  "power- 
off"  situation,  for  reference. 

EXPLOSIVE  CONTAINMENT-CHAMBER.  The  entire  explosive 
electrode  assembly  is  housed  in  a  forty  foot  steel  explosion  chamber 
which  has  been  evacuated  to  1  torr  pressure  for  sound  abatement  of 
the  blast  wave.  Figure  15  is  a  top  view  of  the  entire  experimental 
layout  as  it  presently  exists. 

Diagnostic  Techniques 

OPTICAL  TRACKING  TO  DETERMINE  THE  DETONATION  VELOCITY.  The 
evoluation  of  detonation  enhancement  will  De  followed  with  an  ultra- 
high  speed  image-converter  camera  employing  standard  time- sweep 
techniques.  For  the  sequence  of  experiments  using  the  cylindrical 
electrode  setup,  the  optical  diagnostics  are  performed  while  operating 
in  the  streak  mode  at  a  recording  rate  of  10  mm/ysec.  With  a  possible 
resolution  of  200  nanoseconds,  the  exposure  time  is  10  microseconds 
during  which  the  detonation  front  travels  about  8.8  cm.  The  recording 
sequence  tracks  the  progress  of  the  detonation  front  along  the  length 
of  the  outer  electrode  surface  as  schematically  depicted  in  Figure  16. 
A  Pockel  cell  setup  provides  time-marks. 

The  longitudinal  length  of  the  outer  cylindrical  electrode 
surface  is  vertically  aligned  with  the  primary  focal  plane  of  the 
camera  optics.  The  resultant  image  is  optically  restricted  to  a 
defining  slit  aperature  allowing  only  a  thin  section  of  the  outer 
cylinder  wall  to  be  observed.  Frontal  illumination  is  provided  by 
a  20  microseconds  exploding  wire  source.  This,  in  turn,  is  reflected 
off  the  front-most  section  of  the  outer  circular  cylinder  wall,  a 
portion  of  which  is  subtended  at  the  slit  to  be  re-imaged  and  ampli¬ 
fied  for  recording.  As  the  outer  cylinder  wall  radially  flares,  the 
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FIGURE  13  ELECTRODE  ASSEMBLY:  ANNULAR-COAXIAL  GEOMETRY 
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FIGURE  14  ELECTRICAL-TO-CHEMICAL  feNERGY  COUPLING  IN  THE  CYLINDRICAL 
ELECTRODE  CONFIGURATION. 
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FIGURE  16  STREAK  TRACKING  CONFIGURATION  OF  THE  OUTER  CYLINDRICAL  ELECTRODE, 
SET  UP  TO  FOLLOW  THE  EVOLUTION  OF  DETONATION. 
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reflected  light  is  no  longer  aligned  with  the  slit.  The  subsequent 
"smear"  thus  serves  to  delineate  the  rate  at  which  the  cylinder  wall 
moves  once  the  camera  has  been  calibrated.  The  cylinder  wall  has 
been  optically  polished  and  micro-grooved  every  5.0  mm  down  the 
vertical  extent  of  the  cylinder  wall.  These  precisioned  markings 
serve  to  determine  the  spatial  extent  along  the  run  of  the  cylinder 
wall  to  within  tolerances  of  one- thousandth  of  an  inch.  The  detonation 
velocity  of  PBX  9404  incased  in  the  cylinder  can  be  accurately  measured 
prior  to  an  "enhanced"  run.  The  results  under  "power-on"  conditions 
are  then  photographed  and  measured  from  the  time- sweep  variation  in 
the  luminosity  from  the  reflected  strip  section  off  the  cylinder  wall. 
Subsequently,  the  interrupted  slit  image  of  a  line  in  the  event  is 
swept  across  the  film  emulsion  at  a  uniform  rate  thus  producing  a 
record  of  light  intensity  distribution  against  a  known  time  base. 

This  allows  the  shock  velocity  and  the  duration  to  be  directly  re¬ 
solved  from  a  time-swept  photograph. 

An  alternate  measurement  on  the  coaxial  electrode  is  to 
measure  the  radial  expansion  of  the  cylindrical  electrodes.  Measure¬ 
ments  would  be  made  when  the  electrodes  are  energized  and  when  they 
are  not.  These  two  types  of  measurements  allow  for  relative  compari¬ 
sons  so  enhancement  of  outer  electrode  velocity  can  be  easily  assessed. 

ELECTRICAL  MEASUREMENTS.  Monitoring  the  current  flow 
through  the  electrodes  is  accomplished  by  use  of  a  conventional 
0.2-10  kHz  integrating  Rogowski  Coil  with  a  0.5  mV/A  sensitivity  and 
a  10  nsec,  to-peak  response  characteristic.  The  device  is  designed 
and  shielded  to  be  sufficiently  uncoupled  from  any  external  stray 
noise  sources  such  as,  inductive  pick-up  from  the  storage  coil. 
Calibration  has  been  made  against  a  constant  current  generator  with 
the  same  range  of  signal  characteristics  as  those  expected  to  be 
measur  ed . 


Voltage  measurements  are  made  using  a  fully  calibrated 
resistive  voltage  divider  having  a  characteristic  response  time  on 
the  order  of  about  100  nsec.  The  energy  introduced  into  the  detona¬ 
tion  wave  is  determined  from  a  Tektronix  556  synchronized  scope  trace 
by  the  integration 

W  =  /  T  I{t)  V(t)dt  ,  in  time  =  t. 

*o 

CONDUCTIVITY  MEASUREMENTS.  From  measured  values  of 
the  voltage  drop  across  the  detonation-electrode  assembly  and  of  the 
current  flowing  through  it,  the  total  resistance  can  be  deduced  as 
the  ratio  of  the  two. 29  if  ft  is  tacitly  assumed  that  the  only 
significant  region  of  current  flow  is  through  the  reaction  zone  of 
the  detonation  process,  an  average  value  of  conductivity  can  be 
directly  obtained.  Unfortunately,  conductivity  values  inferred  from 
such  direct  methods  yield  meaningful  results  only  in  special  circum- 


29 

Huddlestone,  R.  H.  and  Leonard,  S.  L. ,  Eds.,  Plasma  Diagnostic 
Techniques ,  Academic  Press,  Inc.,  New  York,  ppT  19-27','  19*35. 
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stances  and  should  be  regarded  with  some  trepidation  when  applied  to 
this  experimental  configuration.  On  the  assumption  that  voltage 
corrections  for  inductive  contributions  from  the  external  circuitry 
can  be  made,  further  complications  due  to  the  possible  existence  of 
a  resistive  layer  caused  by  possible  thermal2^  and  field  emission 
effects^O  in  the  immediate  vicinity  of  the  electrode  are  likely  (even 
under  high  current  conditions).  The  ablative  effects  of  the  detona¬ 
tion  wave  as  it  sweeps  past  the  electrodes  can  produce  a  result 
something  like  a  Chi ld-Langmire  effect.28  The  presence  of  metal 
particles  from  the  electrodes  serve  to  increase  the  conductivity, 
locally  altering  the  electric  field  distribution,  and  consequently 
redistributing  the  potential  across  the  electrode  pair.  Experimental 
evidence  of  these  localized  effects  as  produced  by  condensed  explo¬ 
sive  detonations  are  described  in  the  paper  by  Datsenko,  Korol'kov, 
Kransnov  and  Mel'nikov.22  The  presence  of  the  metal  particles  form 
in  a  layer  near  the  surface  of  the  electrodes. 

Ordinarily,  the  features  of  the  experimental  simplicity  and 
very  localized  sampling  characteristics  would  recommend  electro- static 
probes  as  a  highly  viable  technique  for  conductivity  measurements. 
However,  the  effects  of  a  highly  collision-dominated  medium,  as  in  the 
case  of  the  high  pressure  regime  considered  here,  has  a  two-fold 
consequence.  First,  these  effects  complicate  probe  analysis,  espe¬ 
cially  for  non- equilibrium  situations.  Second,  for  such  high  plasma 
parameter  regemes  (i.e.,  g~12)  the  probe  itself  has  strong  perturba- 
tional  effects  due  to  the  negligible  sheath  formation  under  such 
condi tions.30,31  Finally,  the  latest  pulsed  probe  techniques,  even 
those  employed  in  rapid  fluctuation  measurements,  have  insufficient 
sampling  resolution  (on  the  order  of  ~0.1  ms.). 32  in  conclusion, 
Langmuir  type  probe  techniques  may  not  be  reliable  for  fast  non¬ 
equilibrium  high  pressure  situations  and  the  error  would  be  signifi¬ 
cant  when  some  fraction  of  the  current  discharge  flows  into  the 
probe. 30  The  use  of  various  other  electrostatic  types  of  probes  under 
very  low  current  conditions  have  been  traditionally  employed  in 
attempting  to  measure  conductivity  profiles  in  a  detonation  wave.  In 
these  probe  systems  where  the  instrument  is  mounted  "on  axis"  and 
normal  to  the  detonation  front,  measurements  are  obtained  as  the 
detonation  is  swept  over  them.  These  have  met  with  varying  degrees 
of  success;  the  perturbational  influence  of  the  probe  on  the  detona- 

30 

Huddleston,  R.  H.  and  Leonard,  S.  L.,  Eds.,  Plasma  Diagnostic 

Techniques,  Academic  Press,  Inc.,  New  York,  pp.  113-124,  151-162, 

189-194. 
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Boyd,  R.  F.,  "Probe  Theory  in  the  High  Pressure  Regime,"  Proc.  Phys. 

Soc. ,  (London)  p.  976,  July  1959. 
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Schoenberg,  K.  F.,  "Pulsed  Electrostatic  Probes  as  a  Diagnostic  for 
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tion  wave  itself,  not  to  mention  on  the  conductivity  profile,  and  the 
accorded  resolution  of  such  schemes  are  being  critically  scrutinized 
at  this  time. 

OO 

Microwave  reflection  techniques  do  not  appear  to  suffer 
from  the  major  drawback  of  significantly  perturbing  the  detonation 
process  or  the  conductivity  profile.  But  recently  published  studies 
on  the  latest  known  results  of  this  method  put  the  spatial  resolution 
limits  on  the  order  of  about  3.5  mm.  Thus  the  probe  has  a  time 
resolution  equal  to  the  space  resolution  divided  by  the  detonation 
velocity  of  about  0.4  microseconds.  The  probe  can  readily  be  oriented 
radially  or  axially  to  the  detonation  front.  The  possibility  of 
employing  this  technique  will  have  to  be  further  investigated  in 
future  work.  Application  of  microwave- transmission  techniques  are 
ruled  out  since  electron  densities  higher  than  10^4  cm“3  do  not 
allow  even  150  GHz  microwaves  to  penetrate. 33  Time  Domain  Reflecto- 
metry  is  presently  being  investigated  as  a  possible  resistivity  probe 
technique.  The  study  of  conductivity  profiles  in  partially- ionized 
dense  detonation  systems  warrants  further  investigation. 


A  SUMMARIZING  DISCUSSION  OF  THE  RESULTS. 

The  present  method  for  electrically  coupling  into  the  explosive 
detonation  entails  passing  a  very  high  electric  current  pulse  through 
it.  Specifically,  the  injection  of  170  kA  will  increase  the  conduc¬ 
tion  zone  temperature  -  ultimately  culminating  in  the  overdriving  of 
the  explosive's  detonation  velocity,  and  subsequently,  increasing  the 
detonation  pressure  and  rate  of  energy  release. 

The  preliminary  facility  arrangement  deploys  about  100  kilojoules 
of  electrical  energy  and  is  expected  to  deliver  up  to  about  2  times 
the  explosive's  chemical  release  energy  through  joule  heating,  under 
optimum  power  delivery  conditions.  "Increases  in  detonation  velocity 
are  to  be  measured  by  fast  time-swept,  photographic  techniques.  The 
current  and  voltage  profiles  are  continuously  monitored  by  employing 
standard  diagnostic  instrumentation. 

The  electrically-induced  augmentation  of  the  detonation  wave 
from  a  condensed  explosive  depends  primarily  on  the  explosive's 
electrical  properties.  Briefly,  these  are:  (a)  the  electrical 
conductivity  distribution  behind  the  detonating  explosive’s  shock 
front,  and  (b)  the  breakdown  voltages  of  the  reacted  gases  in  the 
detonation  product  zone,  and  of  the  unreacted  explosive.  The  extent 
of  coupling  this  electrical  energy  into  the  explosive  depends  on  the 
impedance  characterizing  the  LR  network  assembly. 

The  dominant  mechanism  attributed  to  the  energy  enhancement 
effects  are  principally  ohmic  in  nature.  The  responsible  impedance 
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characteristics  consist  of:  (a)  a  pure  resistance,  which  determines 
the  extent  of  joule  heating,  principally  in  the  reaction  zone,  and 
(b)  an  inductive  reactance  which  determines  the  extent  of  inductive 
heating.  For  the  particular  set  of  experimental  parameters  typically 
attributed  to  chemical  explosive  detonation  the  direct  ohmic  heating 
energy  release  mechanism  is  found  to  be  at  least  an  order  of  magni¬ 
tude  smaller  than  the  explosive's  own  chemical  release  energy  effects. 
Furthermore,  for  the  set  of  electrode  configurations  presently  being 
used  inductive  heating  is  insignificant. 

In  conclusion,  this  program  will  afford  the  opportunity  to 
explore  the  feasibility  of  investigating  the  effects  of  augmenting 
the  chemical  energy  of  a  reacting  explosive  with  some  external  energy 
source.  It  will  provide  an  additional  independent  parameter  to 
influence  basic  explosive  properties.  However,  a  number  of  subtasks 
remain  to  be  finished  prior  to  attempting  an  observation  of  enhance¬ 
ment  effects.  Specifically,  a  continued  effort  must  be  made  to 
bring  this  facility  closer  to  a  fully  operational  status.  Once  this 
is  accomplished,  high  current  level  resistance  measurements  must  be 
made.  Finally,  a  demonstration  of  enhancement  effects  will  be 
pursued. 
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APPENDIX  A 

INDUCED  OHMIC  HEATING  EFFECTS 


An  estimate  of  the  "induced"  ohmic  heating  effects  are  readily 
attained  for  the  1.7  x  10^  A  current  flow  in  the  cylindrically 
concentric,  electrode  pair.  The  effect  primarily  depends  on  the 
magnitude  and  rate  of  induced  magnetic  field  diffusion  into  the 
conductive  region  of  the  detonation  process.  These  values  are  based 
on  the  magnetic  Reynold  numberAl  Rm  and  magnetic  diffusion  time  t. 

To  make  the  analysis  more  tractible,  the  conductivity  n  (£ )  will  be 
assumed  constant  across  the  reaction  zone  (i.e.  dn  =  0). 

d t 

Upon^ insertion  of  Equation  (12)  into  Faraday's  Induction 


(A.l) 

.  (A—  2 ) 

The  curl  of  Ampere's  Law  transforms 


Equation, in  R.M.K.S.  units, 


-* 

3  R  ->  ■+ 

ff  =  -  (VxE) 


x  ^nJ 


into 


B 


n 


-+ 
X  J 


(A.3) 


(A.  4) 


Since  there  is  no  curvature, 


V  (V*B)  =  0. 


(A.  5) 


A1 

Sutton,  G.  W.  and  Sherman,  A.,  Engineering  Magnetohydrodynamics, 
McGraw-Hill  Book  Co.,  New  York,  pp.  125-209,  and  pp.  300-304,  1965. 

A2 

Jackson,  J.  D. ,  Classical  Electrodynamics,  2nd,  Ed.,  John  Wiley 
and  Sons,  Inc.f  pp.  471-479,  1975. 
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Relative  to  the  laboratory  frame  of  reference  Equation  (A. 2)  gives 


=  n_  y  2 

9t 


B  +  V  x(DxB) 


(RMKS)  . 


(A. 6) 


The  first  term  on  the  R.  H.  S.  is  the  familiar  field  diffusion  term; 
the  second  term  is  a  convective  term,  important  for  moving  resistive 
fluids.  If  the  field  lines  penetrate  the  conduction  zone  fast 
enough,  1_  =  n  ,  the  induced  field  effects  cause  ohmic  heating  of 
tB  y0£2  -*■ 

the  fluid.  This  energy  is  due  to  the  presence  of  the  B  field;  the 
energy  lost  per  cm2  in  a  time  t  is  (oj2)t.A3  By  Ampere's  law 


V  x  B 


(A.7) 


The  energy  dissipation  is 


(n  j2)  t 


2  ^ 


(A. 8) 


=  2^“  B ,  y Q  =  1.26  x  10  °  Henry/m  .  (A. 9) 

^  o 

A1 

The  magnetic  Reynolds  number  R  (a  comparison  of  the  R.H.S. 
terms  in  Equation  (A6))  will  indicate  now  easily  the  fluid  slips 
through  the  magnetic  field.  Relative  to  the  laboratory  frame,  the 
ratio  of  these  terms  is 


R  =  D  B/& 

m  n/yD  .  B/z2 


U0DZ 


(A. 10) 


For  our  explosive  detonation  situation  where  a  selected  maximum 
resistivity,  n  =  10^  •  m,  and  D  is  the  PBX  9404  detonation 

velocity  max 

D  =  8.8  x  10 3  m/sec 

&  =  1.5  10“3  m  (width  of  idealized  conduction  zone) 


JGriem,  H.  R. ,  and  Loveberg,  R.  H. ,  Eds.,  Methods  of  Experimental 
Physics;  Plasma  Physics,  Vol.  9,  Part  A,  Academic  Press,  Inc., 
New  York,  pp.  403-418,  1970. 
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Rm  = 


„  (1.26  10~6) (8.8xl03) (1.5  10~3) 


'10 


-8 


1.0* 


(A. 11) 


For  Rm  <<  1,  therefore,  some  diffusion  of  the  induced  magnetic 
field  lines  across  the  conductive  zone  of  the  detonation  wave  is 
anticipated.  To  ascertain  how  significant  the  effect  can  become, 
the  most  optimistic  case  would  be  i  =  imax  =  .17  MA.  This,  however, 
assumes  a  rectangular  current  pulse  about  10  microseconds  long, 

where,  in  fact,  the  actual  current  situation  is  i  =  i  exp  (~  t)  ; 

max 

and  consequently  the  magnetic  field  B*Bmax  exp  (^  t)  where,  R  as 

s 

a  function  of  the  conduction  zone  resistivity  n(e),  is  assumed  constant. 
In  the  0.17  MA  situation  for  the  cylindrical  configuration  with 
r1  =  0.7  cm. 

Bmax  =  (  °  ijnax )  =  4.87  Tesla  .  (A.12) 

\  2irr1  / 

From  Equation  (A. 9) ,  the  maximum  energy  dissipation  per  cm  given  by 


( n  j2)  t 


max 


.  4. 

1  yo1max 
2uq  2xr 


_  ^o^max 
8Ti2r2 


(per  unit 
volume) 


(A.  13) 


The  total  energy  available  to  the  5  cm  long  coaxial  electrode  is 


,  *2. 

(n  3  )  xdr 


(A. 14) 


where  the  volume  T  is  a  cylindrical  shell  5  cm  long,  with  radii 
r1/rQ  =  2.0.  Therefore,  using  the  coaxial  electrode  pair 


U 


(1.26  x  10  H  /m)  (1.7  x  10  A)  (0.05  m) 

4(3.141) 


(0.69) 


=  101  joule  .  (A. 15) 

This  energy  release  is  ~  3  x  103  smaller  than  that  of  the  total  chemical 
release  energy  of  the  detonated  charge.  What  is  significant  here, 
however,  is  the  implications  associated  with  Equation  (A. 7) .  Since 
Rm  <<  1,  the  diffusion  term  is  dominant  in  Equation  (A. 6) . 
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From  Equation  (A. 6)  the  magnetic  diffusion  time  is  given  by 

2 

=  y  a£  t  where  £  is  a  length  characteristic  of 
the  spatial  variation  of  B  relative 
to  an  observer  in  the  conduction  zone. 

If  tb  is  defined  as  the  total  time  for  the  field  to  diffuse  through 
the  conduction  zone  as  it  traverses  the  accelerator  electrodes,  the 
power  dissipated  goes  into  joule  heating.  With  the  present  experi¬ 
mental  arrangement,  the  electrode  length  L  =  0.05  m;  taking  a  nominal 
bulk  conductivity  value  a  =  102  mho/m  (cf  Figure  2), 

2 

xB  =  P0ctL  (sec) 

=  (1.25  x  10"6)  (102)  (0.05)2  (A. 16) 

=  0.31  ys. 


The  detonation  " time-of-flight"  through  the  0.05  m  electrodes  is 


td  (=  L/D)  = 


0.05  m 

8. 8x10 J  m/sec 


(A. 17) 


=  5.7  y sec  . 

Consequently,  the  "interaction  parameter” 
tr 

N  e  —  (*Rm) 

tD 

=  0.05  . 


(A. 18) 


N  <  1  implies  that  there  is  sufficient  time  for  the  annihilation 
of  the  magnetic  flux  lines  to  produce  the  ohmic  heating.  Equation 
(A. 6)  implies  that  for  Br  ~.25  MG  (e.g.,  a  multi-turned  solenoid 
about  a  dielectric  shelly  used  as  an  electrode)  the  induction  heating 
energy  would  be  comparable  to  the  chemical  release  energy  of  the 
detonation,  the  explosive  reaction  zone  acting  as  the  secondary  coil 
of  a  transformer ,A^ 


A4 

Hooke,  W. ,  and  Caldirola,  P.,  Eds.,  Symp.  on  Plasma  Heating  and 
Injection,  International  School  of  Plasma  Physics,  Varenna,  Italy, 
pp.  61-72,  Sep  21  -  Oct  1972. 
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